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I.  Summary 

The  second  year  of  an  experimental  study  of  swept  shock  wave  interactions  with 
turbulent  boundary  layers  is  reported.  This  work  relied  entirely  on  non-intrusive,  laser-based 
optical  flow  diagnostics.  Experiments  were  carried  out  to  define  the  flowfield  structure  of 
fin-generated  interactions  over  a  Mach  number  range  from  2.5  to  3.5  using  the  laser  light- 
screen  flow  visualization  technique.  Further  experimental  esulted  in  accurate  skin  friction 
measurements  in  fin-generated  swept  interactions  by  way  of  the  Laser  Interferometer  Skin 
Friction  Meter.  Finally,  techniques  were  perfected  for  the  assessment  of  flowfield  fluctuation 
levels  using  pulsed-laser  holographic  interferometry.  The  results  of  this  research  have  given 
new  insight  into  the  fin-interaction  flowfield  structure,  which  involves  a  supersonic  jet 
produced  by  shock  wave  refraction  and  impinging  on  the  test  surface.  High  skin  friction 
levels  were  measured  in  the  vicinity  of  this  jet  impingement,  which  compare  well  with 
computational  predictions  carried  out  by  others.  Finally,  the  observed  jet  impingement 
structure  is  shown  to  be  similar  to  that  of  leading-edge  shock  wave  impingement,  which  is 
known  to  produce  serious  pressure  and  heating  loads  on  aerodynamic  surfaces  of  high-speed 
flight  vehicles. 


II.  Research  Objectives 

The  research  objectives  for  the  second  year  of  AFOSR  Grant  86-0082,  as  given  in  the 
Statement  of  Work  of  our  proposal  dated  December  1,  1986,  were  as  follows: 

1)  Explore  the  flowfield  structure  of  swept  interactions  due  to  sharp  fins  at  Mach  3  and 
variable  a  (fin  angle  of  attack)  by  way  of  optical  diagnostics. 

2)  For  the  same  flow  conditions,  apply  pulsed-laser  holographic  interferometry  to  assess 
flowfield  fluctuation  levels. 

3)  Obtain  3-D  skin  friction  measurements  in  these  flows  using  the  Laser  Interferometer  Skin 
Friction  Meter. 

Progress  toward  the  achievement  of  these  research  objectives  is  discussed  in  detail  in 
the  following  section. 

III.  Status  of  the  Research  - 
Accomplishments  and  Progress 

1)  Optical  Flowfield  Structure  Results 

The  flowfield  structure  of  shock  wave/turbulent  boundary  layer  interactions  at  Mach 
numbers  2.5,  3,  and  3.5  has  been  investigated  as  proposed  during  the  second  year  of  the 
subject  Grant.  Sharp,  unswept-leading-edge  fins  were  used  to  generate  swept  planar  shock 
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waves  which  interacted  with  the  two-dimensional,  equilibrium  turbulent  boundary  layer  on  an 
adiabatic  flat  plate.  The  shock  wave  strength  was  varied  by  adjustable  fin  angles  of  attack 
up  to  a  maximum  of  18  degrees.  The  flowfield  structure  was  observed  by  the  laser  light- 
screen  technique  (Ref.  1),  in  which  the  beam  from  a  5-watt  Argon-ion  laser  was  expanded  into 
a  thin  sheet  and  passed  through  the  flow  normal  to  the  swept  shock  wave.  Moisture  was 
added  to  the  airstream  of  the  Penn  State  Supersonic  Wind  Tunnel,  so  that  the  cross-section  of 
the  airstream  became  visible  due  to  Mie  scattering  of  the^asfer  sheet.  The  results  were  then 
photographed,  yielding  quantitative  structural  information  over  the  parameter  ranges  stated. 
Since  laser  light-screen  results  are  inherently  of  low  contrast,  digital  image  processing  has 
been  used  for  contrast  enhancement. 

Some  previously-known  features  of  the  flowfield  structure,  such  as  the  separation 
shock,  "lambda-foot,"  and  triple  point,  were  clearly  revealed.  The  evolution  of  these  features 
with  increasing  shock  strength  (ie  increasing  Mach  number  or  fin  angle)  has  been  quantified. 
Weak  interactions  appear  nearly  featureless,  but  distinct  structural  features  appear  and  evolve 
as  the  shock  strength  increases.  From  these  results  we  have  ascertained  the  following: 


*  clear  images  of  detailed  flowfield  structures  can  be  obtained 

*  these  structures  evolve  with  increasing  shock  strength 

*  supersonic  jet  impingement  is  seen  in  strong  interactions 

*  computations  and  surveys  have  not  revealed  these  features 


An  important  discovery  of  this  work  is  the  flowfield  structure  downstream  of  the 
triple-point  in  strongly  separated  interactions,  which  has  not  been  observed  previously.  An 
example  of  such  a  light-screen  image,  with  contrast  enhanced  by  digital  image  processing,  is 
shown  in  Fig.  1.  (Complete  details  of  this  and  related  flow  images  obtained  in  this  research 
are  given  in  Appendix  A,  which  conveys  the  draft  version  of  a  paper  to  be  presented  at  the 
First  National  Fluid  Dynamics  Congress,  Cincinnati  OH,  July  1988.)  A  supersonic  jet  is  seen 
in  Fig.  1  to  emanate  from  the  "lambda-shock"  structure,  curve  downward,  and  impinge  on  the 
test  surface,  leading  to  hign  local  pressure  and  heating.  This  phenomemon  appears  closely 
related  to  Edney-type  (Ref.  2)  leading-edge  shock  interference.  Earlier  flowfield  models,  such 
as  that  of  Zheltovodov  (Ref.  3)  do  not  show  the  flowfield  structure  as  it  is  here  revealed  to 
exist. 


While  both  flowfield  surveys  and  Navier-Stokes  simulations  have  been  carried  out 
previously  for  interactions  similar  to  that  of  Fig.  1,  neither  has  had  the  resolution  to  reveal 
the  critical  downstream  structure.  For  example,  Figs.  2  and  3  show  the  flowfield  structures  of 
strong  fin  interactions  surveyed  by  Kimmel  (Ref.  4)  and  computed  by  Knight  (Ref.  5).  In  the 
former  case,  one  is  forced  to  conclude  that  contours  interpolated  from  flowfield  surveys  are 
simply  not  well  adapted  to  revealing  structural  details  (at  least,  not  at  the  resolution  shown). 
For  the  latter,  the  coarse  mesh  employed  in  the  flowfield  to  economize  on  computer  storage 
causes  structural  features,  even  including  the  main  shock  wave,  to  be  smeared  into  indistinc¬ 
tion.  Discussions  with  D.  Knight  and  C.  Horstman  of  the  AFOSR  3-D  Shock  Wave/Turbulent 
Boundary  Layer  Interaction  Working  Group  reveal  that  there  is  insufficient  memory  in  the 
available  supercomputers  to  allow  resolving  the  flowfield  to  the  required  level. 
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2)  Laser  Interferometer  Skin  Friction  Results 


Measurements  by  Laser  Interferometer  Skin  Friction  Meter  (LISF)  of  the  Mach  3,  a  = 
10°  and  16°  fin  interactions  have  been  successfully  carried  out  in  the  Penn  State  Supersonic 
Wind  Tunnel.  In  order  to  carry  out  these  experiments  it  has  been  necessary  to  reconfigure 
our  LISF  instrument  to  look  directly  down  on  the  flat  plate  test  surface  from  above.  Also,  a 
replacement  test  section  ceiling  was  designed  and  constructed*  to  contain  the  plexiglass  window 
needed  for  optical  observation  from  above.  Having  made  the  necessary  modifications  to  the 
LISF  instrument  and  tunnel  ceiling,  we  ran  kerosene-lampblack  surface  flow  visualization 
traces  of  the  Mach  3,  a  =  10°  and  16°  fin  interactions  to  recheck  the  flow  quality,  which  was 
found  to  be  quite  good. 


Since  LISF  measurements  are  pointwise  and  time-consuming,  some  consideration  was 
given  tc  the  proper  location  of  the  measuring  positions  on  the  flat  plate.  First  of  all,  since 
no  skin  friction  data  of  established  accuracy  exist  for  any  swept  interaction ,  we  wished  to 
produce  data  which  will  serve  as  a  "benchmark.”  It  thus  seems  reasonable  to  acquire  data  in 
the  region  from  the  undisturbed  boundary  layer  to  the  fin  in  the  zone  of  quasiconical  flow 
symmetry,  ie,  outside  the  interaction  inception  zone.  (Measurements  inside  the  inception  zone 
provide  less  general  information,  thus  are  assigned  a  lower  priority.)  Secondly,  we  are  thus 
making  use  of  the  inherent  symmetry  of  the  interaction,  which  calls  for  a  spherical  polar 
system  of  coordinates.  For  surface  measurements,  such  coordinates  form  concentric  circular 
arcs  centered  about  the  virtual  concial  origin  of  the  flow  (see,  eg.  Ref.  6  for  a  detailed 
discussion  of  optimum  interaction  coordinate  systems).  However,  it  was  our  intention  to 
measure  interactions  of  different  a  in  this  way,  and  the  virtual  origin  changes  with  a.  Since 
the  virtual  origin  is  not  far  displaced  from  the  fin  leading  edge,  we  have  chosen  the  latter  as 
the  origin  of  coordinates  for  simplicity.  Note  that  some  investigators  still  remain  confused 
about  the  optimum  coordinate  frame  of  swept  interactions.  We  believe  strongly  that  the 
measurement  of  data  along  lines  parallel  to  the  freestream  in  swept  interactions,  as  was  done 
in  early  studies,  is  confusing  and  inappropriate  in  the  light  of  current  understanding. 


The  chosen  datum  radius  for  detailed  skin  friction  data  in  the  Mach  3,  a  =  10°  fin 
interaction  was  4.5  inches  from  the  fin  leading  edge.  This  was  chosen  to  be  well  downstream 
of  the  inception  zone  according  to  the  established  scaling  law  for  that  zone  (Ref.  1): 


(Lj/50)Re50a  =  *j«tanAs 


where  a  =;  1/3  at  Mach  3,  =;  1133,  and  As  is  the  sweepback  angle  of  the  fin  shock.  For 

our  current  tests  at  Re  =  5.9x1 07/meter  and  50  =  4.2  mm,  Eqn.  1  yeilds  an  inception  length  Lj 
of  about  2  inches  from  the  fin  leading  edge.  By  making  use  of  the  quasiconical  symmetry  of 
the  interaction,  it  is  possible  to  plot  the  LISF  data  vs.  angle  0  from  the  freestream  direction. 
The  results  will  thus  automatically  be  in  dimensionless  coordinates  (as,  indeed,  all  conical 
flows  must  be). 


The  LISF  skin  friction  data  for  the  Mach  3,  a  =  10°  and  16°  interactions  are  shown  in 
these  coordinates  in  Figs.  4  and  5.  Also  shown  are  the  Navier-Stokes  compuational  predictions 
of  Horstman  (private  communication)  and  Knight  (Ref.  7),  Figs.  4  and  5  are  very  revealing 
and  bear  some  discussion. 


First  of  all,  the  reason  to  obtain  LISF  data  is  that  there  is  no  other  known  technique 
which  will  measure  skin  friction  with  known  accuracy  in  a  shock /bound ary  layer  interaction. 
The  Preston -tube  method  is  much  simpler  to  use  than  LISF,  but  does  not  actually  measure  skin 
friction.  Instead,  it  measures  pitot  pressure  at  the  wall  and  infers  skin  friction  through  a 
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calibration  which  cannot  be  trusted  in  a  strongly-interacting  flow.  Further,  reliable  skin 
friction  data  are  sorely  needed  to  test  the  predictive  capabilities  of  current  Navier-Stokes 
solvers.  Here,  representative  state-of-the-art  computations  are  tested  against  reliable  data 
with  an  established  accuracy  of  ±  5%. 

In  brief,  it  appears  that  both  computations  hayj^done  well.  Both  Knight’s  and 
Horstman’s  CFD  results  are  in  good  agreement  with  the-  LISF  skin  friction  measurements  at 
both  fin  angles  of  attack.  Upon  close  inspection,  the  LISF  data  show  a  slight  initial  rise  in 
Cf  at  the  beginning  of  the  interaction,  which  is  displaced  downstream  in  Horstman's  computed 
prediction  and  missing  altogether  from  Knight’s.  Nonetheless,  the  comparison  is  well  sum¬ 
marized  by  a  quote  from  Prof.  Knight  upon  seeing  these  results:  "I’m  surprised  that  the 
computations  did  this  well." 

Though  the  LISF  technique  measures  skin  friction,  in  principle,  with  equal  facility  in 
or  out  of  an  interacting  flow,  we  have  experienced  the  following  problem.  As  shown  in  Figs. 
4  and  5,  Cf  rises  by  as  much  as  a  factor  of  4  in  the  interaction,  above  its  undisturbed  flat- 
plate  level.  At  very  high  wall  shear  levels  the  problem  of  surface-wave  instabilities  on  the  oil 
film  limits  the  number  of  obtainable  LISF  fringes,  and  thus  the  accuracy  of  the  result. 
Monson  (.Ref.  8)  reckoned  the  maximum  wall  shear  level  measureabie  by  LISF  to  be  about  150 
N/m*-.  By  way  of  technique  refinements  and  improved  data  analysis  we  have  reached  perhaps 
four  times  this  value,  though  at  the  cost  of  a  somewhat  larger  error  band  at  peak  Cf  levels. 
A  complete  discussion  of  this  issue,  as  well  as  of  the  LISF  instrument  and  the  experiments 
described  above,  is  included  in  AIAA  Paper  88-0497  which  forms  Appendix  B  of  this  Report. 

In  summary,  we  have  learned  the  following  thus  far  from  the  LISF  experiments: 

*  the  LISF  technique  works  in  swept  interactions 

*  a  surface-wave  limitation  is  present  but  not  overriding 

*  initial  results  compare  well  with  CFD  predictions 

*  more  comparisons  for  stronger  interactions  are  needed 

Finally,  the  high  peak  shear  levels  revealed  in  Figs.  4  and  5  are  connected,  we  believe 
with  the  flowfield  structure  results  described  earlier.  Peak  Cf  is  seen  to  occur  in  the  vicinity 
of  the  supersonic  jet  impingement  noted  in  the  previous  section.  Thus  much  of  the  skin 
friction  behavior  within  these  swept  interactions  may  be  thought  of  as  a  decay  from  peak 
levels  beneath  the  jet  impingement,  forward  to  the  upstream  influence  line.  The  fact  that  the 
computations  predict  this  well  seems  to  indicate  that  the  important  flow  physics  is  present  in 
the  computations,  though  they  do  not  adequate  resolve  such  flow  phenomena  as  the  jet 
impingement. 

Finally,  the  Reynolds-averaged  Navier-Stokes  results  of  Horstman  and  Knight  do  not 
address  flowfield  unsteadiness,  which  seems  almost  certain  to  be  present  in  the  peak  Cf  region. 
For  that  matter,  no  experimental  data  on  unsteadiness  are  available  there  either.  This  is 
definitely  an  issue  for  future  research. 
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3)  Pulsed -Laser  Holographic  Interferometry  Results 

Our  efforts  to  apply  pulsed-laser  holographic  interferometry  to  swept  interactions,  as 
indicated  in  the  earlier  statement  of  research  objectives,  are  still  incomplete  at  the  time  of 
this  writing.  We  have,  however,  completed  the  proof-of-concept  optical  bench  tests  w’hich  are 
prerequisite  to  applying  this  flow  diagnostic  in  actual  wind^nnel  experiments.  A  wind  tunnel 
entry  is  scheduled  for  May,  1988,  to  obtain  holographic  interferograms  in  fin-generated 
shock/boundary  layer  interactions. 

Pulsed-laser  holographic  interferometry  is  both  a  valuable  and  a  difficult  flow  diagnos¬ 
tic.  Its  value  lies  in  its  ability  to  obtain  true  3-D  flowfield  images,  to  generate  a  series  of 
holographic  interferograms  revealing  different  ranges  of  flowfield  fluctuations,  to  image  the 
flow  through  poor  optical  components  or  even  against  opaque  backgrounds,  and  to  provide 
quantitative  density  measurements.  Some  of  the  difficulties  include:  complex  laser  operation, 
limited  lase-  reliability  and  repeatability,  extended  user  training,  and  limited  optical  access  to 
3-D  flows  in  conventional  wind  tunnels. 

We  have  made  significant  progress  in  most  of  these  areas,  with  a  few  difficulties  left 
tc  overcome.  Two  members  of  our  research  team  are  now  proficient  in  the  operation  of  the 
Gas  Dynamics  Lab’s  300  mJ  double-pulse  ruby  laser,  and  others  are  learning  the  technique. 
Beyond  this,  we  have  also  developed  and  demonstrated  the  following  capabilities: 

*  High-quality  holographic  interferogram  production 

*  Both  2-D  and  3-D  interferometric  imaging  of  interacting  flows 

*  Flow  fluctuation  analysis  by  variable-pulse-separation  holography 

*  Flow  unaging  through  plexiglass  and  against  opaque  backgrounds 


Two  examples  of  this  capability  are  shown  here  in  Figs.  6  and  7.  Both  are  photo¬ 
graphs  obtained  recently  in  our  Laboratory,  depicting  reconstructions  of  holographic  inter¬ 
ferograms  of  an  axisymmetric  supersonic  free  jet  (Fig.  6)  and  a  converging  nozzle  with 
internal  shock/boundary  layer  interaction  and  separation  (Fig.  7).  In  both  cases,  small-scale, 
benchtop  flows  were  established  for  technique  development  and  evaluation  (both  flow  pheno¬ 
mena  appear  at  near  full  scale  in  these  Figures).  These  image-plane  holograms  are  recon¬ 
structed  in  white  light  -  one  of  several  possible  methods  -  which  results  in  superimposed 
rainbow  diffraction  colors  on  the  interferograms. 

In  Fig.  6,  an  overexpanded  free  jet  with  Mach  number  about  2.5  exists  from  the 
axisymmetric  nozzle  appearing  in  the  upper  left  of  the  photograph.  This  double-pulse  holo¬ 
graphic  interferogram  was  exposed  once  with  no  flow,  then  again  with  flow  on.  The  entire 
range  of  mean  and  fluctuating  densities  in  the  jet  is  thus  revealed.  Other  interferograms  in 
this  series  (not  shown  here)  were  taken  with  progressively  shorter  intervals  between  laser 
pulses  while  the  flow  was  on.  These  reveal  the  location  and  magnitude  of  flow  fluctuations 
within  a  given  frequency  passband.  For  example,  a  pulse  separation  of  500  fis  results  in  a 
holographic  interferogram  revealing  only  those  flow  components  with  frequencies  greater  than 
1/(500x10-^)  =  2  kHz.  As  this  frequency  high-pass  cutoff  is  raised,  first  the  shock  waves  and 
eventually  the  turbulent  structures  disappear  from  the  interferogram.  For  the  jet  in  Fig.  6, 
pulse  separations  less  than  about  15  /is  reveal  no  remaining  fluctuations  (above  a  67  kHz). 


Fig.  7  is  a  photo  of  the  white-light  interferogram  reconstruction  of  a  "2-D"  supersonic 
nozzle  flow.  Shock-induced  separation  and  an  irregular  shock  crossing  are  apparent  just 
before  the  nozzle  exit  at  the  far  right  in  the  Figure.  The  main  significance  of  this  image  is 
that  it  was  taken  through  two  ordinary  plexiglass  sidewalls.  Internal  strains  in  the  plexiglass 
usually  prevent  other  optical  flow  imaging  methods  (eg  schlieren  and  shadowgraph)  from 
functioning.  We  intend  to  put  this  capability  to  use  by  usy^g.  extended  plexiglass  sidewall  and 
diffuser  windows  in  the  Penn  State  Supersonic  Wind  Tunnel  in  order  to  obtain  the  highlv- 
oblique  viewing  angles  necessary  to  image  swept  shock/boundary  layer  interactions.  These 
windows  are  already  fabricated  and  are  awaiting  use  in  an  upcoming  tunnel  entry. 

The  remaining  problems  to  be  solved  in  implementing  holographic  interferometry  as 
described  in  the  Research  Objectives  section  mainly  concern  the  pulsed-ruby  laser  itself. 
Although  it  is  essentially  new,  it  suffers  occasional  operational  lapses  which  have  delayed  its 
use  in  actual  wind  tunnel  tests  until  now.  We  are  currently  negotiating  with  the  instrument’s 
manufacturer  (Apollo  Lasers,  Inc.)  to  resolve  these  problems. 

Meanwhile,  as  just  discussed,  we  have  developed  an  arsenal  of  holographic  inter¬ 
ferometry  techniques  and  capabilities  during  the  grant  year  reported  here  that  will  allow  us  to 
accumpiisn  ihis  part  of  the  stated  research  objectives.  We  regard  this  worx  as  a  top  priority 
item,  insofar  as  time-resolved  data  in  swept  interactions  are  extremely  sparse,  and  none  of  the 
available  results  address  the  fluctuations  of  the  interaction  flowfield  at  all.  Thus  the  time- 
resolved  nature  of  these  flow  is  almost  completely  unknown.  It  is,  further,  a  key  issue  in  our 
fundamental  understanding  of  swept  interaction  behavior,  as  well  as  in  the  validation  of 
current  CFD  predictions  which  are  time-averaged.  Experimental  results  directly  addressing 
this  issue  are  anticipated  during  the  first  quarter  of  the  1988-89  AFOSR  grant  year. 


4)  Other  Significant  Accomplishments 

In  addition  to  the  research  described  above,  carried  out  in  direct  fulfillment  of  the 
Statement  of  Work  for  the  second  year  of  the  subject  Grant,  additional  related  work  has  been 
done  during  the  second-year  reporting  period.  In  particular,  we  have  put  into  operation  the 
research  equipment  obtained  under  AFOSR  Grant  84-184,  a  DoD-URIP  Grant  now  expired  which 
was  specifically  for  equipment  in  support  of  our  research  on  shock/boundary-layer  and  related 
high-speed  viscous-inviscid  interactions. 

An  example  is  documented  in  Appendix  C,  which  contains  AIAA  Paper  88-2062  on  the 
subject  of  microprocessor  control  of  high-speed  wind  tunnel  stagnation  pressure.  The 
microprocessor  controller  and  transducers  obtained  under  AFOSR  Grant  84-184  were  put  into 
operation  to  provide  automatic  stagnation  pressure  control  of  the  Penn  State  Supersonic  Wind 
Tunnel  Facility.  The  accompanying  analysis  and  experimental  verification  have  also  resulted  in 
a  Masters  Degree  for  one  member  of  our  research  group. 

Many  other  examples  could  be  cited  of  equipment  obtained  earlier  under  AFOSR  Grant 
84-184  which  is  now  in  use  and  contributing  to  the  currently-reported  research.  These  items 
include,  in  particular,  the  test  section  assembly.  Argon-ion  laser,  Chronolite-8  and  Hycam  II 
high-speed  cameras,  and  data  acquisition  system  obtained  under  the  DoD-URIP  program.  These 
are  cited  here  because  of  their  crucial  impact  on  our  ability  to  have  obtained  the  successful 
results  described  in  this  report. 


IV.  Plans  for  Third-Year  Effort 


The  third  year  of  AFOSR  Grant  86-0082  is  now  under  way.  Our  goals  for  this  effort 
center  around  the  use  of  our  established  capability  in  non-intrusive,  laser-based  optical  flow 
diagnostics  to  attempt  to  resolve  some  of  the  key  questionscpncerning  the  swept  shock/boun¬ 
dary  layer  interaction  phenomenon.  Early  in  this  year.-^’D  pulsed-laser  holographic  inter¬ 
ferometry,  as  described  above,  will  be  brought  to  bear  on  the  problem  of  swept  interaction 
unsteadiness  and  structure.  Our  further  goals  are  listed  below,  from  the  Statement  of  Work 
of  our  proposal  dated  October,  1987: 

1)  Further  investigate  the  detailed  structures  of  swept  fin-generated  interactions,  includ¬ 
ing  both  weak  and  strong  cases,  using  pulsed-laser  holography. 

2)  Apply  laser  light-screen  diagnostics  to  examine  the  growth  of  weak  and  strong 
interactions  with  distance  downstream  of  the  interaction  apex. 


3)  Widen  the  current  database  of  LISF  skin  friction  measurements  in  these  flows  to 
include  Mach  numbers  higher  than  3. 

4)  Explore  the  improvement  of  LISF  instrument  performance  at  high  wall  shear  levels. 


V.  Publications 


1)  Settles,  G.  S.,  Metwaily,  O.  M.,  Hsu,  J.,  and  McIntyre,  S.  S.,  "Visualization  of  High-Speed 
Flows  at  the  Penn  State  Gas  Dynamics  Laboratory,"  Invited  paper  to  be  presented  at  the  1988 
International  Conference  on  Lasers  and  Electro-Optics,  Santa  Clara,  Oct.  30,  1988.  (in 
preparation.) 

2)  Fung,  Y-T,  Settles,  G.  S.,  and  Ray,  A.,  "Microprocessor  Control  of  High-Speed  Wind  Tunnel 
Stagnation  Pressure,"  fo  be  presented  at  the  AIAA  Aerodynamic  Testing  Conference,  San 
Diego,  CA,  May  18-20,  1988. 

3)  Lu,  F.  K.,  and  Settles,  G.  S.,  "Structure  of  Fin-Shock/Boundary  Layer  Interactions  by  Laser 
Light-Screen  Visualization,"  to  be  presented  at  the  1st  National  Fluid  Dynamics  Congress, 
Cincinnati,  July  1988. 

4)  Kim,  K-S,  and  Settles,  G.  S.,  "Skin  Friction  Measurements  by  Laser  Interferometry  in  Swept 
Shock  Wave/Turbulent  Boundary- Layer  Interactions,"  AIAA  Paper  38-0497,  presented  at  the 
AIAA  26th  Aerospace  Sciences  Conference,  Reno,  NV,  January  11-14,  1988.  ( Submitted  for 
publication  to  AIAA  Journal  and  now  undergoing  review.) 

5)  Settles,  G.  S.,  "Aerospace  and  Wind  Tunnel  Testing,"  Chapter  IV-4  of  Handbook  of  Flow 
Visualization ,  ed.  W.-J.  Yang,  Hemisphere  Press.  (To  be  published  in  1988.) 

6)  Kim,  K.-S.,  and  Settles,  G.  S.,  "Skin  Friction  Measurements  by  Laser  Interferometry,"  Ch.  3 
of  an  AGARDograph  entitled  A  Survey  of  Measurements  and  Measuring  Techniques  in  Rapidly 
Distorted  Compressible  Turbulent  Boundary  Layers,  eds.  H.  H.  Fernholz,  A.  J.  Smits,  and  J.-P. 
Dussauge,  (To  be  published  in  1988.) 
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7)  Lu,  F.  K.,  Settles,  G.  S.,  and  Horstman,  C.  C.,  "Mach  Number  Effects  on  Conical  Surface 
Feat.,  ss  of  Swept  Shock  Boundary-Layer  Interactions,"  AIA A  Paper  87-1365,  June  1987, 
<  I  '  jepted  for  publication  in  AIAA  Journal.) 


VI.  Personnel 

Principal  Investigator 

Prof.  Gary  S.  Settles 


Graduate  Research  Assistants 

Joseph  C.  Hsu  (Ph.D.  candidate  supported  by  the  subject  Grant) 

Kwang-Soo  Kim  (Ph.D.  candidate  supported  by  the  subject  Grant.  Dissertation  entitled  "The 
Laser  Interferometer  Skin  Friction  Meter  for  Compressible  Flows,"  in  preparation. 
Graduation  expected  during  1988.) 

Frank  K.  Lu  (Ph.D.  candidate  supported  in  part  by  the  subject  Grant  and  in  part  by  NASA- 
Ames  Interchange  Agreement  NCA2-192.  Thesis  title:  "Fin-Generated  Shock  Wave/- 
Boundarv  Layer  Interactions."  Degree  will  be  granted  on  May  14,  1988.) 


VII.  Interactions 

1)  Spoken  Papers  at  Meetings,  Conferences,  and  Seminars 


1)  Kim,  K-S,  and  Settles,  G.  S.,  "Skin  Friction  Measurements  by  Laser  Interferometry  in  Swept 
Shock  Wave/Turbulent  Boundary-Layer  Interactions,"  AIAA  Paper  88-0497,  presented  at  the 
AIAA  26th  Aerospace  Sciences  Conference,  Reno,  NV,  January  11-14,  1988. 

2)  Settles,  G.  S.,  "Comments  on  Useful  vs.  Unproductive  Code  Validation  Exercises,"  presented 
at  the  NASA  CFD  Validation  Workshop,  Moffett  Field,  CA,  July  14-16,  1987. 

3)  Settles,  G.  S.,  "Current  and  Planned  Shock  Wave/Turbulent  Boundary  Layer  Research  at 
Penn  State,"  presented  at  the  First  Annual  3-D  Shock/Boundary  Layer  Interaction  Working 
Group  Meeting,  University  Park,  PA,  July  27-28,  1987. 

4)  Lu,  F.  K.,  Settles,  G.  S.,  and  Horstman,  C.  C.,  "Mach  Number  Effects  on  Conical  Surface 
Features  of  Swept  Shock  Boundary-Layer  Interactions,"  AIAA  Paper  87-1365,  presented  at  the 
AIAA  19th  Fluid  and  Plasma  Dynamics  Conference,  Honolulu,  HI,  June  1987. 
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2)  Consultative  and  Advisory  Functions 


•  Dr.  C.  C.  Horstman  of  NASA-Ames  Research  Center  visited  our  group  on  April  6-7,  1987. 
Detailed  discussions  cf  shock/boundary  layer  interaction  problems  were  held. 

•  Dr.  A.  G.  Havener  of  SUNY-Binghamton  University  v2i^5d  our  group  on  April  9,  1987.  His 
expert  advice  on  holographic  interferometry  issues  was  sought. 

•  PI  visited  Princeton  University  May  7,  198"7,  to  participate  as  External  Reader  in  the  Ph.D. 
Defense  of  Roger  L.  Kimmel,  a  student  of  Prof.  Bogdonoff  whose  work  concerned  swept 
interactions. 

•  PI  and  F.  K.  Lu  attended  AIAA  19th  Fluid  and  Plasma  Dynamics  Conference,  Honolulu,  HI, 
June  1987.  AFOSR-supported  work  was  presented  (see  above).  Many  discussions  were  held 
with  colleagues  involved  in  swept  interaction  research  and  related  issues. 

•  PI  and  several  group  members  attended  AFOSR/ONR  Contractors  Meeting  at  Penn  State  the 
week  of  June  21,  1987.  Discussions  were  held  with  Dr.  Tishkoff  of  AFOSR  and  many  others. 

•  PI  visited  US  Naval  Weapons  Center,  China  Lake,  CA,  on  July  13,  1987  to  assist  Drs.  E. 
Gutmark  and  K..  Schadow  with  optical  flow  diagnostic  experiments. 

•  PI  attended  NASA  CFD  Validation  Workshop  at  Ames  Research  Center  July  14-16,  1987,  and 
presented  AFOSR-supported  research  to  participants  (see  above). 

•  PI  participated  in  transonic  wing  tests  at  NASA-Ames  Research  Center,  July  14,  1987,  and 
instructed  personnel  there  on  the  proper  use  of  surface  flow  visualization  method. 

•  PI  met  with  R.  S.  Boray  of  AFWAL/POPT  July  16,  1987  and  discussed  LISF  skin  friction 
instrument.  Copy  of  technical  paper  was  later  sent  to  Boray. 

•  PI  met  with  Lockheed  Corp.  personnel  T.  Harsha,  L.  Long,  and  P.  Bevilacqua  on  July  17, 
1987  to  discussed  cooperation  on  hypersonic  swept  interaction  experiments. 

•  PI  discussed  with  F.  Staub  of  GE  Corp.  (Schenectady,  NY)  possible  PSU  facility  use  for 
govt.-sponsored  GE  testing  program,  July  23,  1987. 

•  PI  hosted  the  First  Annual  3-D  Shock/Boundary  Layer  Interaction  Working  Group  Meeting, 
University  Park,  PA,  July  27-28,  1987.  Composed  primarily  of  AFOSR-supported  investigators, 
this  Group  produced  a  document  outlining  future  research  needs.  Dr.  J.  D.  Wilson  of  AFOSR 
was  present  along  with  8  others. 

•  Dr.  E.  Gutmark  of  the  US  Naval  Weapons  Center,  China  Lake,  CA,  visited  our  group  on  Aug. 
12,  1987,  presented  a  seminar,  and  participated  in  detailed  research  discussions. 

•  PI  discussed  with  E.  Kawecki  of  Pratt  &  Whitney  Corp.  possible  consulting  arrangement  for 
advice  to  P&W  NASP  program  on  shock/boundary  layer  interactions,  Sept.  25,  1987. 

•  Dr.  M.  Reichman  of  ONR  visited  our  group  on  Sept.  30,  1987.  He  was  given  a  tour  of  our 
facilities  and  discussions  of  high-speed  propulsion  research  were  held. 
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•  Mr.  George  Donahue  of  the  RAND  Corp.  -risked  Penn  State  on  Sept.  30,  1987.  PI  discussed 
with  him  issues  concerning  NASP  aerodynamics  and  propulsion. 

•  Dr.  C.  C.  Horstman  of  NASA-Ames  Research  Center  visited  our  group  on  October  26-27, 

1987.  Detailed  discussions  of  shock/boundary  layer  interaction  problems  were  held. 

•  Dr.  S.  Lekoudis  of  ONR  visited  Penn  State  on  October-^?,*  1987  regarding  proposed  super¬ 
sonic  mixing  research  under  joint  ONR/AFOSR  sponsorship. 

•  Dr.  T.  Doligalski  of  the  US  Army  ARO  visited  our  group  on  October  28,  1987.  Discussions 
were  held  of  mutual  research  interests  and  a  facility  tour  was  given. 

•  Dr.  A.  Gross  of  NASA-Ames  visited  Penn  State  on  December  3,  1987,  and  was  given  a  tour 
of  our  facilities  and  a  presentation  on  our  research  program. 

•  PI  was  quoted  in  NY  Times  article  January  5,  1988,  concerning  swept  interaction  research 
and  high-speed  test  facilities. 

•  PI  and  K-S  Kim  attended  AIAA  26th  Aerospace  Sciences  Meeting,  Reno,  NV,  Jan.  I !  - 1 4 , 

1988.  Presented  AFOSR-supported  work  (see  above).  Attended  sessions  on  shock/boundary 
layer  interactions  and  related  topics.  Held  discussions  with  co-workers  in  this  research  area 
including  S.  M.  Bogdonoff,  C.  C.  Horstman,  D.  D.  Knight,  D.  S.  Dolling,  &  G.  R.  Inger  at  2nd 
Meeting  of  the  3-D  Shock/Boundary  Layer  Interaction  Working  Group. 

•  PI  met  with  Dr.  J.  McMichae!  of  AFOSR  to  discuss  current  research  program  and  plans,  Jan. 
12,  1988  in  Reno,  NV. 

•  PI  received  invitation  to  attend  Soviet  Conference  on  Simulation  Problems  and  Wind  Tunnels 
from  Dr.  A.  A.  Zheltovodov  (ITAS  Novosibirsk)  on  Jan.  19,  1988  (Declined  due  to  conflict  with 
1st  Natl.  Fluid  Dynamics  Congress.) 

•  PI  discussed  heat  transfer  instrumentation  with  R.  Neumann  of  AFWAL  on  Feb.  24,  1988. 

•  PI  discussed  possible  visiting  scientist  arrangement  at  AFWAL  with  J.  Sliang  and  K.  Stetson 
of  AFWAL  on  March  8,  1988  (at  their  request;  negotiations  in  progress). 

•  PI  rendered  advice  to  Rocketdyne  Corp.  personnel  J.  Skipper,  J.  Blauer,  and  L.  Rubin  on 
holographic  interferometry  for  NASP  testing  program,  March  10,  22,  and  24,  1988. 
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Pressure  Contours  of  the  M,*,  =  3,  a  =*  20'  Fin  Interaction,  from  Flowfield 
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Structure  of  Fin-Shock/Boundary-Layer  Interactions 
by  Laser  Light- Screen  Visualization 

F.  K.  Lu*  and  G.  S.  Settlesf,  Pennsylvania  State  University,  University  Park,  PA,  16802 


Abstract 

The  flowfield  structure  of  fin- generated  shock- wave/turbulent  boundary-layer  interac¬ 
tions  at  nominal  Mach  numbers  of  2.5,  3,  and  3.5  is  studied  using  laser  light-screen  flow 
visualization.  A  sharp-leading-edged  fin  was  used  to  generate  swept  planar  shock  waves 
which  interacted  with  the  2-D,  equilibrium  turbulent  boundary  layer  on  an  adiabatic  flat 
plate.  The  shock  wave  strength  was  varied  by  adjustable  fin  angles  up  to  18° .  The  beam  of 
a  4  W  Argon-ion  laser  was  expanded  into  a  thin  sheet  and  passed  through  the  flow  roughly 
normal  to  the  interaction  sweep  line.  Moisture  in  the  airstream  allowed  the  photography  of 
interaction  “slices.”  Quantitative  data  on  the  flowfield  structures  were  obtained  at  Mach 
3  and  3.5.  Weak  interactions  appear  nearly  featureless,  but  distinct  structural  features 
appear  and  evolve  as  the  shock  strength  is  increased.  An  important  discovery  of  this  work 
is  the  flowfield  structure  downstream  of  the  shock  triple-point,  which  has  not  been  prop¬ 
erly  observed  before.  A  supersonic  jet  bounded  by  the  slip  line  and  flow  separation  zone  is 
seen  to  curve  downward  and  impinge  on  the  flat  plate.  This  phenomenon  appears  closely 
related  fc  Edney-type  leading-edge  shock  interference.  Such  structural  features  have  not 
been  resolved  in  previous  flowfield  surveys  or  Navier-Stokes  solutions. 


Nomenclature 

Mrj0  =  incoming  freestream  Mach  number 

pQ  =  stagnation  pressure  of  incoming  stream,  MPa  (psia) 

Re  =  unit  Reynolds  number,  m~ 1  (/ft) 

Taw  =  adiabatic  wall  temperature,  K  (°  R) 

Tw  =  wall  temperature,  K  (°  R) 

T0  =  freestream  stagnation  temperature,  K  (°R) 

a  =  angle  made  by  fin  with  respect  to  the  incoming  freestream  direction,  deg. 

(3  =  angle  made  by  surface-flow  features  or  flowfield  features  projected  onto  the  test 

surface,  with  respect  to  the  incoming  freestream  direction,  deg. 

0O  =  angle  made  by  inviscid  shock-wave  trace  on  the  test  surface  with  respect  to  the 
incoming  freestream  direction,  deg. 

8  =  local,  undisturbed  boundary-layer  thickness,  mm  (in.) 

Subscripts 

A1  =  primary  flow  attachment  line  (Figure  10) 

B  =  top  of  separated  boundary  layer  between  S2  and  S3  (Figure  6b) 
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e  =  edge  of  vortical  structure  (Figure  6b) 
h  =  height  of  flowfield  feature  above  test  surface,  mm  (in.) 

I  =  distance  of  flowfield  feature  projected  onto  the  test  surface  from  the  fin  apex 
p  =  intersection  of  the  separation  shock  with  the  boundary  layer  (Figure  6b) 

52  =  separation  shock  (Figure  6b) 

53  =  internal  shock  (Figure  6b) 

t  =  top  of  vortical  structure  (Figure  6b) 

tp,  TP=  triple  point  of  lambda  shock  structure  (Figure  6b) 

U  =  upstream  influence  line 

<t>  —  angle  mads  by  straightline  segment  of  flowfield  structure  in  the  plane  normal  to 

the  inviscid  shock 

oo  =  incoming  freestream  conditions 

Introduction 

The  problems  of  shock  wave /boundary- layer  interactions  have  been  prominent  in  mod¬ 
ern  fluid  dynamics  for  several  decades.  These  problems  not  only  pace  the  development  of 
advanced  computational  and  theoretical  methods,  but  also  have  important  relevance  to 
high-speed  flight.  Numerous  publications  and  several  survey  articles  have  appeared  on  the 
subject  (see,  e.  g.  ,  the  recent  review  of  swept  interactions  by  Settles  and  Dolling1). 

Most  recently,  shock  wave/boundary-layer  interaction  problems  have  taken  on  added 
importance  due  to  proposals  for  the  development  of  a  new  generation  of  supersonic  and 
hypersonic  vehicles.  Although  these  vehicles  will  likely  be  designed  by  computational 
methods,  the  lack  of  previous  experience  in  computational  hypersonics  reduces  confidence 
in  its  application.2  This  is  especially  so  in  the  computation  of  complicated  viscous-inviscid 
interacting  fiowfields  such  as  that  of  a  shock-wave/boundary-layer  interaction3  . 

Thus,  to  support  computational  code  validation  and,  more  importantly,  to  understand 
the  fluid-dynamic  phenomena  involved,  recourse  to  experiments  on  shock/boundary-layer 
interactions  is  necessary.  To  meet  the  first  goal,  a  broad  and  detailed  database  is  called 
for  in  order  to  verify  evolutionary  efforts  in  code  development.  The  second  goal  requires 
detailed  experiments  as  well,  but  not  always  the  same  experiments.  Exploratory  studies 
aimed  at  revealing  the  interaction  behavior,  extent,  and  structure  may  be  more  useful  than 
detailed  quantitative  measurements  in  satisfying  this  goal.  Such  experiments  often  should 
cover  wide  parametric  ranges,  say,  of  Mach  and  Reynolds  numbers. 

The  present  paper  reports  such  an  experimental  study,  aimed  at  understanding  the 
structure  of  the  interaction  of  a  fin-generated,  swept  shock  wave  with  a  turbulent  bound¬ 
ary  layer.  This  basic  research  configuration,  which  is  depicted  in  Figure  1,  represents 
practically-occurring  interactions  on  high-speed  vehicles.  Two  examples  are  fin-fuselage 
junctions  and  corner  flows  in  air  inlets. 

A  key  traditional  question  in  all  three-dimensional  shock/boundary-layer  interactions 
is  the  flowfield  structure.  For  the  often-studied  fin  case  of  Figure  1,  previous  experiments 
including  detailed  yaw-probe  surveys4-7  have  provided  some  information  regarding  this 
structure,  as  sketched  in  Figure  2  (this  sketch,  from  Reference  1,  depicts  the  interaction 
flowfield  in  a  plane  normal  to  the  inviscid  shock).  In  interactions  of  sufficient  strength, 
the  incoming  boundary  layer  is  separated  by  a  “separation  shock”  which  is  oblique  to  the 
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oncoming  flow  component  in  the  normal  plane.  This  shock  is  therefore  relatively  weak 
■  ^mpared  to  the  main,  freestream  shock.  Hence  another,  downstream,  shock  is  needed  to 
turn  the  interacting  flow  in  order  to  satisfy  the  boundary  conditions  imposed  by  the  fin. 
This  system  of  shock  waves  is  similar  to  that  seen  in  the  more  familiar  case  of  the  two- 
dimensional  transonic  shock/boundary  layer  interactk^i..  When  visualized  using  optical 
techniques,  such  a  shock  system  appears  to  have  a  “A”  shape,  hence  the  term  lambda-foot 
structure  is  used  to  describe  it.  At  the  “triple  point”  of  the  lambda-foot,  where  the  three 
shocks  converge,  a  slip-surface  forms.  Across  this  surface,  which  appears  as  a  line  in  the 
plane  of  Figure  2,  pressure  and  flow  direction  must  match  but  all  other  flow  quantities  may 
be  discontinuous.  Beneath  the  lambda-foot,  the  separated  boundary  layer  rolls  up  into  a 
flattened  vortex  as  indicated  in  Figure  2.  In  some  cases  a  secondary-separation  vortex  is 
suspected,  although  there  is  still  insufficient  direct  evidence  to  support  this. 

While  these  flowfield  features  are  relatively  well-understood,  what  happens  down¬ 
stream  of  the  lambda-foot  is  not  so  clear.  Figure  2  shows  an  expansion  and  recompression 
based  roughly  on  the  inferences  of  previous  investigators.4,7  However,  the  phenomenol¬ 
ogy  of  this  interaction  region  has  not  been  clearly  shown,  either  by  previous  flowfield 
surveys1,4-6  or  by  computational  simulations  of  the  flowfield.8  In  all  cases,  insufficient 
resolution  was  obtained  in  previous  studies  to  clearly  portray  the  downstream  structure. 
In  the  present  work  we  will  demonstrate  by  flow  visualization  that  the  earlier  model  of  the 
downstream  compression-expansion  region,  as  shown  in  Figure  2,  is  incorrect. 

Flowfield  visualization  is  not  new  to  shock/boundary-layer  interactions.  It  has  been 
applied  in  several  previous  studies.9- 14  Of  these,  References  9  through  11  used  light-screen 
visualization,  Reference  12  used  the  electron-beam  technique,  and  References  13  and  14 
used  conical  shadowgraphy.  In  general,  the  three-dimensional  nature  of  swept  interactions 
poses  difficulty  in  applying  classical  optical  techniques.  However,  for  dimensionless  fin- 
generated  interactions,  quasiconical  symmetry15  allows  some  simplification  of  both  the 
experimental  technique  and  the  data  analysis.  Further,  while  most  of  the  cited  studies 
sought  only  qualitative  flow  visualization  results,  a  proper  analysis  could  have  been  used 
to  extract  quantitative  information  from  the  images  obtained.16  Finally,  the  previous 
applications  of  light-screen  flow  visualization  to  fin-generated  shock-wave/boundary-layer 
interactions9- 11  shared  the  common  problem  of  low  contrast  and  resolution,  making  it 
very  difficult  to  extract  quantitative  information  from  them. 

Hence,  to  fulfill  the  need  for  a  better  understanding  of  the  interaction  flowfield  struc¬ 
ture,  a  new  set  of  laser  light-screen  observations  of  fin-generated  interactions  has  been 
performed  and  is  reported  here.  The  aims  of  this  study  are  to  examine  the  flowfield  struc¬ 
ture  in  detail,  and  to  extract  quantitative  data  on  the  key  structural  features  from  these 
observations.  Before  describing  these  experiments  it  is  appropriate  to  briefly  review  the 
laser  light-screen  technique  which  they  employed. 


The  Loser  Light-screen  Technique 


Light-screen  (also  light-sheet  and  vapor-screen)  flow  visualization  is  a  standard  exper¬ 
imental  technique  originating  from  the  early  1960’s.17  Though  early  difficulties  were  expe¬ 
rienced  in  collimating  non-coherent  light,  the  recent  availability  of  powerful  laser  sources 
now  enables  high-quality  images  to  be  obtained  without  undue  difficulty.  Usually  the 
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flow  is  seeded  with  water  vapor,  either  naturally  or  deliberately  by  injecting  water  into  the 
airstream18,  though  other  seeding  materials  have  been  tried  as  well.  Light  is  Mie-scattered 
by  the  aerosol  particles  as  the  flow  passes  through  the  light-screen. 

Different  scattering  intensities  thus  indicate  different  flow  densities,  assuming  that 
the  aerosol  was  initially  evenly-distributed  in  the  flowj^However,  boundary-layers,  if  in 
contact  with  relatively  hot  surfaces,  may  re-evaporate  "condensing  moisture,  making  inter¬ 
pretation  of  the  results  at  such  locations  difficult.  Further,  a  high-speed  boundary-layer, 
being  less  dense  than  the  freestream,  contains  a  smaller  particle  number-density  than  the 
freestream.  It  is  thus  observed  that  such  boundary-layers  appear  dark  in  light-screen  im¬ 
ages.  Also,  boundary-layer  fluid  which  leaves  the  surface  during  flow  separation  appears 
dark  as  well.17, 18 

One  drawback  of  the  light-screen  technique  arises  from  condensing  moisture  in  a  su¬ 
personic  stream,  which  adversely  affects  the  flow.19  For  example,  condensation  shocks  in 
the  expansion  of  a  supersonic  wind  tunnel  nozzle  are  known  to  alter  the  flow  properties, 
sometimes  seriously.  Many  previous  studies  have  used  moist  air  with  a  frostpoint  as  high 
as  — 4  3  C  (25  3F),  which  jeopardizes  the  results  due  Lo  the  possibility  of  condensation  shock 
formation.19  The  procedure  we  have  followed  to  prevent  this  from  happening  is  described 
in  the  next  section. 

The  light-screen  technique  has  been  popular  for  visualizing  the  leeside  flows  over  high¬ 
speed  missiles  and  delta  wings  at  angle-of-attack.  For  further  details  on  the  technique  and 
such  applications  thereof  in  high-speed  flows,  the  reader  is  referred  to  McGregor’s  original 
work17  and  the  more  recent  contribution  of  Snow  and  Morris18. 

The  Experiment 

Wind  Tunnel  and  Test  Model 

The  Penn  State  Supersonic  Wind  Tunnel,  used  in  this  study,  is  a  blowdown  facility 
with  a  nominal  Mach  number  range  of  1.5  to  4.  The  variable  Mach  number  capability  is 
achieved  by  way  of  an  asymmetric  sliding-block  nozzle.  The  test  section  of  the  wind  tunnel 
is  150  mm  (6  in.)  wide,  165  mm  (6.5  in.)  high  and  610  mm  (24  in.)  long.21,22 

For  the  present  experiments  the  interaction  test  surface  (Figure  1)  is  a  flat  plate, 
500  mm  (19.5  in.)  long,  which  spanned  the  test  section.  A  fin  model  with  a  10°  sharp 
leading  edge  was  placed  with  its  tip  216  mm  (8.5  in.)  from  the  plate  leading  edge  and 
26.2  mm  (1.03  in.)  from  the  tunnel  sidewall.  The  fin  was  100  mm  (4  in.)  high,  127  mm 
(5  in.)  long  and  6.35  mm  (0.25  in.)  thick.  The  fin  height  of  about  30<5  ensured  that 
the  interaction  was  a  semi-infinite  one.23  The  length  of  the  fin  was  chosen  to  provide 
the  maximum  interaction  extent  while  allowing  sufficiently  large  angles-of-attack  to  be 
obtained  without  stalling  the  wind  tunnel.  Further,  the  height- to- length  ratio  ensured 
that  disturbances  from  the  top  of  the  fin  leading  edge  passed  behind  the  test  region. 

The  fin  was  held  tightly  onto  the  flat  plate  by  a  fin-injection  mechanism  mounted 
on  the  tunnel  sidewall.  The  bottom  of  the  fin  had  a  rubber  seal  which  ensured  that  no 
leakage  under  the  fin  occurred  during  the  tests.  Once  test  conditions  were  established,  the 
fin-injection  mechanism  pneumatically  injected  the  fin  to  a  preset  angle-of-attack,  a.  This 
was  necessary  only  for  tests  with  a  larger  than  approximately  14°.  At  lower  angles,  a  was 
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fixed  before  the  run.  The  fin  angle  was  determined  to  0.1°  accuracy  using  a  machinist’s 
protractor. 

Test  Conditions 

The  experiments  were  performed  with  freestream^Iach  numbers  of  2.47,  2.93,  and 
3.44.  The  initial  conditions  of  the  incoming  freestream  at  these  Mach  numbers  are  sum¬ 
marized  in  Table  1. 

The  undisturbed  flat-plate  boundary-layers  at  these  Mach  numbers  have  been  sur¬ 
veyed  and  documented22 .  These  boundary-layers  are  turbulent  and  in  equilibrium  in  the 
sense  that  they  satisfy  the  combined  iaw-of-the-wail/law-of-the-wake  with  constant  wake- 
strength  parameter  within  the  accepted  equilibrium  range.  The  test  boundary- layers  are 
also  near-adiabatic  ( Tw  /Taw  ss  1.04).  Example  velocity  profiles  across  the  undisturbed 
boundary-layers  at  the  two  higher  Mach  numbers  are  shown  in  Figures  3a  and  3b,  along 
with  the  corresponding  Sun-Childs24  wall-wake  curvefits  to  the  data.  The  surveys  along  the 
fiac-place  cencerilne  and  38  mm  (1.5  in.)  to  each  sine  showed  that  the  test  ooundary-layers 
are  two-dimensional. 

The  dewpoint  (or,  more  correctly,  the  frost  point)  of  the  air  in  the  storage  tank  for 
the  present  experiments  was  kept  to  — 14  ±  2°C  (6  ±  5°F),  measured  when  the  storage- 
tank  pressure  was  1.4  MPa  (200  psia).  This  dewpoint  results  in  an  error  due  to  moisture 
condensation  of  less  than  one  percent  of  the  freestream  static  pressure,  or  the  freestream 
Mach  number,  for  the  Mach  number  range  of  the  present  experiments.21-22  This  amount 
of  water  vapor  was  also  sufficient  to  yield  good  light-screen  images.  In  order  to  achieve 
this  moisture  level,  the  air,  after  being  compressed,  is  sent  directly  to  the  storage  tank 
without  passing  through  the  normal  drying  equipment.  Experience  with  the  Penn  State 
Supersonic  Wind  Tunnel  reveals  that  the  residual  amount  of  water  vapor  in  the  air  after 
compression  results  in  a  dewpoint  of  about  — 14°  C. 

Discrete  fin  angles-of-attack  were  set  at  a  =  7° ,  10° ,  15°  and  18° .  However,  at 
Moo  =  2.47,  the  a  =  18°  case  was  found  to  cause  tunnel  stalling  and  was  thus  eliminated 
from  the  test  matrix. 

Optical  Equipment 

The  experimental  optical  setup  for  laser  light-screen  imaging  is  diagrammed  in  Fig¬ 
ure  4.  The  collimated  beam  of  a  4  W  Argon-ion  laser  is  reflected  from  two  steering  mirrors, 
SMl  and  SM 2,  before  passing  through  a  cylindrical  lens  which  spreads  the  beam  into  a 
sheet  of  light  1.2  mm  (0.05  in.)  thick.  The  steering  mirrors  allow  for  quick  and  easy  adjust¬ 
ment  uf  the  beam.  By  rotating  mirror  SMl  and  translating  and  rotating  mirror  SM 2,  the 
light-screen  could  be  positioned  accurately  within  limits  imposed  by  the  300  mm  (6  in.) 
mirror  diameters. 

The  light-screen  was  photographed  by  an  Olympus  OM-2S  camera  with  a  Vivitar 
28-90  mm  macrofocusing  zoom  lens  combined  with  a  Vivitar  2x  macrofocusing  telecon¬ 
verter.  Photographs  were  taken  automatically  and  were  bracketed  by  ±2/3  and  ±4/3 
stops,  yielding  five  photographs  per  test  run.  Using  Kodak  TMAX  ISO-400  film,  the 
optimum  exposure  time  was  found  to  be  in  the  range  of  1/15  to  1/30  second. 
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To  provide  a  reference,  a  “tare”  photograph  was  taken  of  a  grid  of  2.54  mm  (0.1  in.) 
squares.  This  grid  was  positioned  to  coincide  with  the  light-screen  plane  in  the  wind  tunnel 
test  section.  During  subsequent  data  reduction  the  geometry  of  the  Sow  structures  was 
obtained,  through  triangulation,  by  referring  to  the  tare  photograph. 

The  Sn  and  Sat  plate  were  painted  Sat-black  to  reduce  stray  reSections.  Nonetheless, 
with  the  light-screen  incident  on  the  Sn,  reSected  light  caused  photographs  from  the  initial 
experiments  to  have  poor  contrast.  This  situation  was  remedied  by  passing  the  light-screen 
about  5  mm  (0.25  in)  to  the  rear  of  the  Sn.  Most  of  the  incident  laser  light,  then  reSected 
from  the  far  wall  of  the  wind  tunnel,  no  longer  entered  the  camera  lens  (Figure  5).  This 
orientation  of  the  light-screen,  however,  prevented  it  from  being  exactly  normal  to  the 
undisturbed  Sn  shock  wave,  so  that  a  further  correction  was  required  in  the  data  reduction. 
Finally,  the  camera  was  positioned  approximately  normal  to  the  light-screen  in  order  to 
capture  as  much  side-scattered  light  as  possible. 

The  results  of  the  experiments,  a  set  of  photographic  negatives,  were  printed  on  high- 
contrast  photographic  paper.  A  selection  of  these  results  was  then  digitized  using  a  high- 
resolution  video  camera  feeding  a  Quantex  QX-9210  digital  image  processor.  Contrast 
enhancement  and  psuedocolor  encoding  were  performed  on  the  480  x  640-pixel  digitized 
images.  Though  not  attempted  thus  far,  it  appears  possible  to  correlate  the  digitized  image 
brightness  level  with  the  flow  density,  utilizing  the  known  flow  density  levels  before  and 
after  the  fin  shock  as  reference  conditions. 

Results 

Qualitative  Analysis 

Based  on  Korkegi’s  well-known  separation  criterion,25  incipient  three-dimensional  tur¬ 
bulent  boundary-layer  separation  occurs  at  Mach  2.5,  3,  and  3.5  for  shock  waves  generated 
by  a  =  7°,  6°,  and  5°  fins,  respect ively.  Therefore  me  present  results  are  generally  ex¬ 
pected  to  involve  boundary-layer  separation,  which  is  likely  to  be  quite  strong  for  the 
larger  fin  angles. l5~9i3~  15  21-23  One  thus  expects  to  observe  flowfield  features  similar  to 
those  sketched  in  Figure  2.  To  illustrate  this  point,  Figure  6a  presents  a  digitally-enhanced 
photograph  of  the  light-screen  visualization  obtained  at  Mach  3.5  with  an  a  —  10°  fin.  The 
essential  features  are  labelled  in  the  schematic  sketch  shown  in  Figure  6b. 

The  photograph  of  Figure  6a  reveals  that  the  freestream  shock  ST  is  curved  just 
above  the  triple-point,  TP.  (Wide-field  photographs  show  th^t  51  is  eventually  straight 
further  away  from  its  interaction  with  the  boundary-layer.)  The  incoming  boundary  layer, 
BL,  separates  beneath  the  initial  oblique  separation  shock,  52,  according  to  the  accepted 
flowfield  model.  The  separation  shock  compresses  and  slows  the  streamtube  subtended 
by  the  points  SP  and  TP.  The  internal  or  rear  shock,  53,  is  also  obvious  in  the  light- 
screen  image.  Thus,  the  well-known  lambda-foot  structure  of  shock-induced  boundary- 
layer  separation,  as  described  previously,  can  be  clearly  seen  in  the  present  results  typified 
by  Figure  6a.  However,  no  clear  structure  is  visible  in  the  photograph  downstream  of  the 
lambda-foot  (that  is,  to  the  left  of  it  in  Figure  6a).  Our  contention,  elaborated  below,  is 
that  a  definite  structure  is  present  but  too  weak  in  this  example  to  be  visible. 

Shock  waves  52  and  53  are  essentially  straight,  as  is  the  “top”  of  the  boundary  layer 
between  SP  and  ST.  Thus  the  angles  <^S2,,  <j>s3  and  <f>B ,  in  the  plane  normal  to  the 


inviscid  freestream  shock,  can  be  defined  for  these  features  and  can  be  extracted  from  the 
light-screen  photographs  in  the  form  of  quantitative  data. 

At  M-n  =  3.5  and  with  a  increased  to  15° ,  flowfield  structural  features  downstream  of 
the  lambda-foot  are  now  strong  enough  to  be  visible  in  the  laser  light-screen  image  (Figure 
7a).  In  particular,  a  curved  line  (presumably  the  slip-l'uje^emanates  from  the  triple-point 
and  impinges  on  the  flat  plate.  A  high  brightness  (density)  level  is  obvious  between  the  slip 
line  and  the  dark  boundary-layer  separation  region.  Further,  striations  are  visible  normal 
to  the  curved  boundaries  of  this  downstream  flow  region.  These  features  are  also  visible  in 
photographs  taken  at  =  3.5  and  a  =  18° .  Prior  evidence  of  the  structure  revealed  here 
in  the  interaction  flowfield  downstream  of  the  lambda-foot  is  known  to  the  authors  from 
Mach  6  interaction  light-screen  photographs  by  T.  J.  Goldberg  of  NASA-Langley  Research 
Center.  Goldberg’s  work  remains  unpublished,  but  was  sketched  and  briefly  discussed  by 
Oskam,  et.  al.11 

Figure  7a  is  highly  reminiscent  of  the  leading-edge  shock  interference  phenomenon 
pioneered  by  Edney.26  It  resembles  a  “Type  IV”  shock  interference  Edney’s  parlance,  which 
is  an  interference  resulting  in  a  supersonic  jet  impingement  upon  the  adjacent  surface  (see 
Figure  8).  By  analogy  with  Edney’s  results,  this  jet-impingement  phenomenon  is  believed 
to  shed  new  light  on  the  high  surface  heating  rates  which  have  been  observed  near  flow 
reattachment  in  separated,  swept  shock/boundary-layer  interactions.27 

Specifically,  based  on  the  above  observations,  we  propose  a  revised  model  for  the 
interaction  structure  downstream  of  the  lambda-foot.  Instead  of  the  previous  model  shown 
in  Figure  2,  that  of  Figure  7b  now  appears  more  reasonable.  The  streamtube  that  is 
processed  through  lambda-shocks  52  and  53  is  bounded  downstream  by  the  slip-line  “e” 
(above)  and  the  vortical  separation  region  “h”  (below).  This  streamtube  is  forced  to  bend 
downward  toward  the  test  surface  by  the  high  pressure  in  region  “n”  after  the  inviscid  shock 
“c.”  Prandtl-Meyer  expansion  fans  “f”  accelerate  the  streamtube  toward  the  surface  in  the 
form  of  a  supersonic  jet.  Finally,  a  normal-shock  “g”  occurs  just  above  the  impingement 
of  the  jet  upon  the  test  surface  in  the  vicinity  of  the  flow  attachment  line  “1.” 

The  new  interaction  structure  model  of  Figure  7b  is  presently  supported  only  for 
strong  interactions  with  significant  three-dimensional  boundary-layer  separation.  For  those 
weaker  interactions  in  the  present  test  matrix,  the  photographic  light-screen  results  do  not 
clearly  show  any  downstream  structure.  However,  direct  visual  observations  of  the  light- 
screen  during  the  experiments  did  reveal  the  slip-line  shape  shown  in  Figure  7b  in  the  much 
weaker  interaction  at  =  3  and  a  =  10°.  This  is  in  direct  contrast  to  the  previous 
interpretation  by  Oskam,  et.  al.4  of  his  flowfield  survey  results  from  the  same  interaction. 
(Oskam  concluded  that  the  boundary-layer  did  not  separate.)  We  believe  that  the  available 
flowfield  survey  data  are  of  insufficient  resolution  to  reveal  the  features  shown  in  Figures 
7a  and  7b,  or,  for  that  matter,  to  reveal  boundary-layer  separation.  For  the  record,  the 
available  Navier-Stokes  solutions  for  these  interactions,  typified  in  Reference  8,  similarly 
fail  to  resolve  the  flowfield  structure  because  of  insufficient  grid  resolution. 

Finally,  Figure  9  shows  an  additional  light-screen  image  for  the  weak  interaction  case 
of  =  2.5  and  a  =  10° .  Here,  the  curvature  of  the  fin  shock  is  still  seen  but  its 
bifurcation  into  a  lambda-foot  is  quite  weak.  From  an  examination  of  Figures  6a,  7a,  and 
9,  one  observes  that  the  rear  leg  of  the  lambda-foot  rotates  downstream  (clockwise)  as  the 


interaction  strength  increases,  thus  decreasing  the  angle  <f>S3.  It  further  appears  from  this 
comparison  that  the  interaction  flowfield  structure  is  a  strong  function  of  the  interaction 
strength,  that  is,  the  degree  of  boundary-layer  separation.  These  qualitative  observations 
are  pursued  further  in  the  discussion  which  follows. 


Quantitative  Analysis 

Previous  studies18-13-15,21-23  have  shown  that  the  swept- inter  act  ion  flowfield,  away 
from  an  inception  region  near  the  fin  tip,  appears  to  grow  conically  from  a  virtual  origin. 
This  principle  was  dubbed  quasiconical  interaction  symmetry  by  its  originators.15  Figure  10 
illustrates  this  concept  in  terms  of  features  observed  in  earlier  surface-flow  visualization21 
experiments.  For  the  present  experiments,  quasiconical  symmetry  is  exploited  in  the  data 
analysis.  The  points  SE ,  ST,  SP  and  TP  of  Figure  5b  are  projected  onto  the  flat  plate 
from  the  light-screen  photographs,  whence  the  angles  of  these  features  with  respect  to 
the  flow  direction  can  be  measured.  (The  angular  reference  scheme  is  diagrammed  in 
Figure  11.1  For  simplicity,  these  measurements  are  referenced  to  the  fin  leading-edge 
rarher  than  the  true  conical-flow  virtual  origin,  which  introduces  a  certain  error.  This 
error  is  believed  to  be  generally  comparable  to  the  angular  accuracy  limit  of  about  ±1.5° 
due  to  the  resolution  of  the  light-screen  photographs.  In  some  cases,  however,  the  leading- 
edge-origin  simplification  may  lead  to  more  serious  errors. 

The  separation  shock  S'?,  is  used  to  illustrate  how  d>s;.  <j>S3  and  d>B  are  determined. 
52  intersects  the  flat  plate  at  C  in  the  plane  of  the  lightsheet  Pi.  In  the  plane  normal 
to  the  inviscid  shock,  P 2,  52  intersects  the  flat  plate  at  C' ,  yielding  LBC' A  =  d>S2.  Let 
A3  =  htp  and  AO  =  ltp.  The  line  OCC  makes  an  angle  e  with  respect  to  the  incoming 
freestream.  Note  that  Pi  and  P2  intersect  along  AB.  Therefore,  from  trigonometry, 


<ps:  =  arctan 


|_LP  tan  (e  -  0tp) 


In  general,  any  straight-line  segment  can  be  reso’-'-ed  in  a  similar  manner,  so  that 


<p  =  arctan 


l  tan  |e  —  0\ 


The  absolute  value  |e  —  0\  is  used  here  so  that  <j>  <  90° . 

Angular  Characterization  of  Interaction  Structure 

Having  used  the  preceding  analysis  to  extract  quantitative  information  from  the  laser 
light-screen  results  at  Mach  2.95  and  3.44,  we  may  now  examine  certain  aspects  of  the 
interaction  structure.  Note  that,  according  to  the  principle  of  quasiconical  interaction 
symmetry,  it  makes  no  sense  to  characterize  the  structure  in  terms  of  either  spatial  positions 
or  length  dimensions.  Only  angles  are  necessary  or  desirable  for  such  a  description. 

The  angle  made  by  the  downstream  edge  of  the  separated,  vortical  region  with  respect 
to  the  incoming  freestream  direction,  /?,,  is  plotted  against  the  fin  angle,  a,  in  Figure  12. 
Also  plotted  for  comparison  in  Figure  12  is  the  attachment-line  angle,  /?A1,  obtained  from 
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previous  surface-flow  visualization  results.21,22  It  can  be  seen  that  is  slightly  smaller 
than  /?A1  at  both  Mach  numbers,  though  these  two  angles  exhibit  similar  trends  with  a. 
The  error  incurred  earlier  in  approximating  the  virtual  conical  origin  by  the  fin  leading-edge 
position  in  determining  0e  is  believed  to  be  responsible  for  this  discrepancy.  Physically, 
0e  >  /?a  i  would  be  expected. 

Next,  the  intersection  of  the  separation  shock  withUfre  boundary  layer  makes  an  angle 
f3p  with  the  freestream,  when  projected  onto  the  test  surface  (see  Figure  11)  The  angle 
j3p  is  plotted  in  Figure  13  against  a.  Also  shown  for  comparison  in  the  Figure  is  the 
upstream-influence  angle,  ,  reported  previously.21  It  can  be  seen  from  this  plot  that, 
within  the  data  scatter, 

/ %  —  P u 

Thus,  one  may  conclude  that  the  upstream- influence  line  coincides  with  the  origin  of  the 
compression  system  which  becomes  the  separation  shock-wave.  This  makes  physical  sense, 
and  has  been  observed  often  in  two-dimensional  interactions  by  comparing  schlieren  photos 
and  wall  pressure  distributions. 

The  triple-point  of  the  lambda-foot  shock  structure  is  usually  assumed  to  be  coplanar 
with  the  undisturbed,  inviscid  shock  wave.  However,  in  view  of  the  shock  curvature  dis¬ 
cussed  earlier,  this  cannot  be  strictly  true.  Figure  14  reveals  this  by  comparing  the  growth 
of  the  projected  triple-point  angle  on  the  test  surface,  f3tp ,  with  the  inviscid-shock  angle 
P0.  At  both  Mach  2.93  and  3.44  the  triple-point  is  from  0.5  to  2°  ahead  of  the  inviscid 
shock.  Further,  Figure  14  indicates  that  the  triple-point  angle  approaches  the  inviscid 
shock  angle  with  increasing  shock  strength  (meaning  that  the  inviscid  shock  curvature  is 
reduced).  The  physical  explanation  of  this  is  not  immediately  obvious. 

Finally,  the  separation-shock  angle  <j)s i  is  plotted  versus  fin  angle  a  in  Figure  15. 
A  dramatic  decrease  in  d>S2  is  seen  as  the  interaction  grows  stronger.  Physically,  this 
corresponds  to  a  stretching  of  the  interaction  “footprint”  over  a  wider  angular  spread. 
Assuming  locally  two-dimensional  flow  in  the  plane  normal  to  the  inviscid  shock,  <j>S2  is 
specified  if  <j>B  and  the  normal  component  of  freestream  Mach  number  are  known.  There 
is  previous  evidence  from  free-interaction  arguments  and  supporting  experiments  that  the 
flow  separation  angle  («  d>B )  remains  relatively  fixed  between  about  8°  and  12°  in  many 
two-  and  three-dimensional  supersonic  interactions.1  Thus,  the  observed  decrease  of  (f> S2 
in  Figure  15  appears  to  be  due  primarily  to  the  increase  of  the  normal  component  of 
freestream  Mach  number  with  increasing  a. 

Conclusions 

In  summary,  this  flow-visualization  study  has  examined  the  structure  of  fin-generated 
shock- wave/ turbulent  boundary- layer  interactions  over  ranges  of  both  Mach  number  and 
interaction  strength.  The  resulting  laser  light-screen  images  have  permitted  a  clear  ob¬ 
servation  of  the  interaction  structure  as  well  as  the  extraction  of  quantitative  data  on 
structural  features.  Several  observations  and  conclusions  have  been  drawn  from  this  study. 
Principally,  it  has  been  shown  that  a  previously  unknown  jet-impingement  structure  exists 
downstream  of  the  familiar  lambda-shock  of  well-separated  interactions.  This  suggests 
the  phenomenological  reason  for  the  high  heat  transfer  rates  observed  by  others  in  the 
flow  reattachment  region,  in  direct  analogy  with  the  problem  of  leading-edge  shock  wave 


interference.  By  contrast,  however,  weaker  interactions  without  significant  flow  separation 
appear  relatively  featureless  except  for  a  noticeable  curvature  of  the  main  shock  wave.  Pre¬ 
vious  interaction  studies  involving  both  flowfield  surveys  and  Navier-Stokes  computations 
have  failed  to  resolve  the  structural  features  observed  here  by  optical  means. 

Acknowledgments*^* 

This  research  was  supported  by  AFOSR  Grant  86-0082  from  the  U.  S.  Air  Force  Office 
of  Scientific  Research,  monitored  by  Dr.  J.  D.  Wilson,  and  by  Joint  Research  Interchange 
NCA2-192  with  the  NASA  Ames  Research  Center,  in  which  Dr.  C.  C.  Horstman  was  the 
research  collaborator. 

References 

1.  Settles,  G.  S.  and  Dolling,  D.  S.,  “Swept  Shock  Wave/Boundary-Layer  Interactions,” 
in  AIAA  Progress  in  Astronautics  and  Aeronautics:  Tactical  Missile  Aerodynamics, 
edited  by  M.  Hemsch  and  J.  Nielsen,  Vol.  104,  AIAA,  New  York,  1986,  pp.  297-279. 

2.  Anderson,  J.  D.,  Jr.,  “A  Survey  of  Modern  Research  in  Hypersonic  Aerodynamics,” 
AIAA  Paper  84-1578,  June  1984. 

3.  Horstman,  C.  C.,  “Prediction  of  Hypersonic  Shock- Wave/Turbulent-Boundary-Layer 
Interaction  Flows,”  AIAA  Paper  87-1367,  June,  1987. 

4.  Oskam,  B.,  “Three-Dimensional  Flow  Fields  Generated  by  the  Interaction  of  a  Swept 
Shock  Wave  with  a  Turbulent  Boundary  Layer,”  Ph.D.  Thesis,  Princeton  University, 
September  1976. 

5.  Peake,  D.  J.,  “The  Three-Dimensional  Interaction  of  a  Swept  Shock  Wave  with  a  Tur¬ 
bulent  Boundary  Layer  and  the  Effects  of  Air  Injection  on  Separation,”  Ph.D.  Thesis, 
Carleton  University,  Canada,  March  1975. 

6.  Kimmel,  R.  L.,  “An  Experimental  Investigation  of  Quasi-Conical  Shock  Wave/- 
Turbulent  Boundary  Layer  Interactions,”  Ph.D.  Thesis,  Princeton  University,  June 
1987. 

7.  Zheltovodov,  A.  A.,  “Regimes  and  Properties  of  Three-Dimensional  Separation  Flows 
Initiated  by  Skewed  Compression  Shocks,”  Zhurnal  Prikladnoi  Mekhaniki  i  Tekhnich- 
eskoi  Fiziki,  No.  3,  May-June  1982,  pp.  116-123. 

8.  Knight,  D.  D.,  Horstman,  C.  C.,  and  Shapey,  B.  L.  “The  Flowfield  Structure  of  the 
3-D  Shock  Wave-Boundary  Layer  Interaction  Generated  by  a  20  Degree  Sharp  Fin  at 
Mach  3,”  AIAA  Paper  86-0343,  January,  1986. 

9.  Kubota,  H.,  “Investigations  of  Three-Dimensional  Shock  Wave  Boundary  Layer  Inter¬ 
actions,”  Ph.D.  Thesis,  Cranfield  Institute  of  Technology,  January  1980. 

10.  Rao,  D.  S.,  Private  Communication  in  Peake,  D.  J.  and  Tobak,  M.,  “Three- 
Dimensional  Interactions  and  Vortical  Flows  with  Emphasis  on  High  Speeds,”  NASA 
TM  81169,  March  1980. 


m 


11.  Goldberg,  T.  J.,  Private  Communication  in  Oskam,  B.,  Vas,  I.  E.  and  Bogdonoff, 
S.  M.,  “Mach  3  Oblique  Shock  Wave/ Turbulent  Boundary  Layer  Interactions  in  Three 
Dimensions,”  AIAA  Paper  76-0336,  January  1976. 

12.  Watson,  R.  D.,  “Experimental  Study  of  Sharp-  and  Blunt-Nose  Streamwise  Corners 
at  Mach  20,”  NASA  TN  D-7398,  April  1974.  - 

13.  Zubin,  M.  A.  and  Ostapenko,  N.  A.,  “Structure  of  Flow  in  the  Separation  Region 
Resulting  from  Interaction  of  a  Normal  Shock  Wave  with  a  Boundary  Layer  in  a 
Corner,”  Izvestiya  Akademii  Nauk  SSSR,  Mekhanika  Zhidkosti  i  Gaza,  Vol.  14,  May- 
June  1979,  pp.  51-58. 

14.  Settles,  G.  S.  and  Teng,  H.-Y.,  “Flow  Visualization  Methods  for  Separated  Three- 
Dimensional  Shock  Wave/Turbulent  Boundary  Layer  Interactions,”  AIAA  Journal, 
Vol.  21,  March  1983,  pp.  390-397. 

15.  Settles,  G.  S.  and  Lu,  F.  K.,  “Conical  Similarity  of  Shock/  Boundary-Layer  Interactions 
Generated  by  Swept  and  Unswept  Fins,”  AIAA  Journal,  Vol.  23,  No.  7,  pp.  1021-1027, 
July  1985. 

16.  O’Hare,  J.  E.,  “Digital  Image  Analysis  for  Aerodynamic  Testing,”  in  Proceedings  of 
the  27th  Symposium  of  ISA,  Vol.  18,  Part  One,  1981. 

17.  McGregor,  I.,  ‘Development  of  the  Vapour  Screen  Method  of  Flow  Visualization  in 
a  3  ft  x3  ft  Supersonic  Tunnel,”  in  AGARDograph  70,  Flow  Visualization  in  Wind 
Tunnels  Using  Indicators,  compiled  by  R.  L.  Maltby,  April  1962. 

18.  Snow,  W.  L.  and  Morris,  O.  A.,  “Investigation  of  Light  Source  and  Scattering  Medium 
Related  to  Vapor-Screen  Flow  Visualization  in  a  Supersonic  Wind  Tunnei,”  NASA  TM 
86290, 1984. 

19.  Wegener,  P.  P.  and  Mack,  L.  M.,  “Condensation  in  Supersonic  and  Hypersonic  Wind 
Tunnels,”  Advances  in  Applied  Mechanics,  Vol.  V,  edited  by  H.  L.  Dryden  and  Th.  von 
Karman,  Academic  Press,  New  York,  1958,  pp.  307-447. 

20.  Monaghan,  R.  J.,  “Tests  of  Humidity  Effects  on  Flow  in  a  Wind  Tunnel  at  Mach 
Numbers  between  2.48  and  4,”  British  A.R.C.  CP  No.  247,  January  1955. 

21.  Lu,  F.  K.,  Settles,  G.  S.  and  Horstman,  C.  C.,  “Mach  Number  Effects  on  Conical 
Surface  Features  of  Swept  Shock  Boundary-Layer  Interactions,”  AIAA  Paper  87-1365, 
June  1987. 

22.  Lu,  F.  K.,  “Mach  Number  Effects  on  Fin-Generated  Shock-Wave  Boundary-Layer 
Interactions,”  Ph.D.  Thesis,  The  Pennsylvania  State  University,  to  be  published. 

23.  McClure,  W.  B.,  “An  Experimental  Study  into  the  Scaling  of  the  Interaction  of 
an  Unswept  Sharp  Fin-Generated  Shock/Turbulent  Boundary  Layer  Interaction,” 
M.S.E.  Thesis,  Princeton  University,  January  1983. 

24.  Sun,  C.  C.  and  Childs,  M.  E.,  “A  Modified  Wall-Wake  Velocity  Profile  for  Turbu¬ 
lent  Compressible  Boundary  Layers,”  AIAA  Journal  of  Aircraft,  Vol.  10,  June  1973, 
pp.  318-383. 


25.  Korkegi,  R.  H.,  “A  Simple  Correlation  for  Incipient  Turbulent  Boundary-Layer  Separa¬ 
tion  due  to  a  Skewed  Shock  Wave,”  AIAA  Journal,  Vol.  11,  November  1973,  pp.  1578- 
1579. 

26.  Edney,  B.  E.,  “Anomalous  Heat  Transfer  and  Pressure  Distributions  on  Blunt  Bodies 
at  Hypersonic  Speeds  in  the  Presence  of  an  Impinging  Shock,”  FFA  Report  115,  The 
Aeronautical  Research  Institute  of  Sweden,  Stockholm,  Sweden,  February  1968. 

27.  Neumann,  R.  D.,  and  Hayes,  J.  R.,  “Introduction  to  Aerodynamic  Heating  Analysis 
of  Supersonic  Missiles,”  in  AIAA  Progress  in  Astronautics  and  Aeronautics:  Tactical 
Missile  Aerodynamics,  edited  by  M.  Hemsch  and  J.  Nielsen,  Vol.  104,  AIAA,  New 
York,  1986,  pp.  421-481. 


Table  1  Incoming  Freestream  Conditions 


Po,  MPa 
(psia) 

T0,  K 

(°R) 

Re  x  10  6 ,  m~ 1 
(/ft) 

2.47  —  0.1% 

0.54x2.0% 

(78) 

295  ±  0.9% 
(531) 

53.8  ±0.9% 

(16.3) 

2.95  ±0.3% 

0.76  ±2.7% 

(110) 

295  ±  0.9% 
(531) 

58.9  ±  1.9% 

(17.8) 

3.44  ±  0.2% 

1.03  ±  3.0% 

(150) 

295  x  0.8% 
(531) 

64.0  ±  1.7% 

(19.4) 

Figure  1 
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Comparison  of  Separation  Shock  Origin  &  Upstream  Influence 
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Abstract 

Experimental  measurements  have  been  made  of  wall  shear  stresses  in  swept  interactions  of  a  planar  shock 
wave  generated  by  a  fin  at  angle  of  attack,  and  the  two-dimensional  turbulent  boundary  layer  on  a  flat  plate. 
Test  tonditicus  were:  Mach  number  3.03,  Reynolds  number  Rtv  =  i4913,  wall  temperature  near  aaiabatic. 
and  fin  angles  10°  and  16°.  Measurements  were  made  using  the  Laser  Interferometer  Skin  Friction  (LISF) 
Meter,  which  optically  detects  the  rate  of  thinning  of  an  oil  film  applied  to  the  test  surface.  The  results 
show  that  such  measurements  are  practical  in  high-speed  interacting  flows,  and  that  a  repeatability  of  ±Q% 
or  better  is  obtainable.  Further,  with  proper  data  acquisition  and  reduction,  LISF  measurements  appear 
feasible  at  very  high  wall  shear  levels  which  were  previously  considered  unobtainable.  Dramatic  increases  in 
wall  shear  were  observed  in  both  swept  interactions  tested.  Data  obtained  along  conical  rays  agree  with  the 
quasiconical  interaction  model  proposed  previously.  The  data  are  compared  with  the  computational  Navier- 
Stokes  solutions  of  Horstman  and  Knight,  wherein  k-e  and  algebraic  eddy-viscosity  turbulence  models  were 
used,  respectively.  Both  computations  simulate  some  features  of  the  measured  :>  distributions  but  fail  to 
predict  other  features  properly. 
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Nomenclature 

=  skin  friction  coefficient  based  incoming  freestream  conditions 
=  incoming  freestream  Mach  number 
=  number  of  fringes 

=  stagnation  pressure  of  incoming  stream,  MPa  (psia) 

=  wail  static  pressure  on  flat  plate,  MPa  (psiaj 
=  incoming  freestream  static  pressure,  MPa  (psia) 

=  radial  distance  measured  from  the  fin  leading-edge,  cm  (in.) 

=  Reynolds  number  based  on  the  local,  undisturbed  boundary-layer  momentum  thickness 

=  distance  measured  along  surface  streamline,  mm 

=  adiabatic  wall  temperature,  K 

=  wall  temperatu  K 

=  freestream  stagnation  temperature,  K 

=  dimensionless  wall-wake  velocity  and  height  coordinates 

=  orthogonal  streamwise,  normal,  and  spanwise  coordinates,  cm 

=  angle  made  by  fin  with  respect  to  the  incoming  freestream  direction,  deg. 

=  angle  with  respect  to  freestream  direction,  measured  from  fin  leading-edge,  deg. 

=  local,  undisturbed  boundary-layer  thickness,  mm 

=  distance  from  oil-film  leading  edge  to  laser  spot  along  surface  streamline,  mm 
=  oil  viscosity,  centistokes 
=  boundary-layer  wake-strength  parameter 
-  wall  shear  stress,  N/m2 


Introduction 

Swept  shock  wave/boundary-layer  interactions  are  among  the  chief  aerothermodynamic  problems  which 
currently  limit  high-speed  propulsion  and  flight.  They  also  represent  a  fundamental  problem  area  of  modern 
fluid  dynamics  which  has,  thus  far,  evaded  any  comprehensive  theoretical  or  computational  treatment.  Many 
publications  and  a  recent  survey  article1  have  appeared  on  this  subject. 

The  computation  of  complex  vbcous-inviscid  interacting  flows  such  as  the  swept  shock/boundary-layer 
interaction  is  presently  paced  by  the  need  for  detailed,  “benchmark*  experiments  for  code  validation.2 
In  particular,  it  is  widely  recognized  that  traditional  measurements  such  as  surface  pressure  and  oil-flow 
visualization  are  inadequate  for  this  purpose.  The  current  state-of-the-art  Navier-Stokes  solvers  still  employ 
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highly-simplified  turbulence  models,  the  efficacy  of  which  is  in  considerable  doubt  for  complex  interacting 
flows.  Experimental  data  which  directly  address  the  turbulent  stresses  in  these  flows  are  sorely  needed  to 
advance  both  the  development  and  validation  of  turbulence  models.  Unfortunately,  such  experiments  are 
extremely  difficult  and  expensive  to  perform  in  high-speed  flows. 

Recognizing  this,  the  participants  of  the  1980-81  AFOSR-HTTM-Stanford  Conference  on  Com¬ 
plex  Turbulent  Flows  scrutinized  the  available  database  and  recommended  certain  directions  for  future 
experimentation.3  Among  these  was  the  suggestion  that  detailed  aftl-accurate  heat  transfer  and  skin  friction 
data  be  obtained  in  high-speed  flow  experiments.  These  quantities,  being  directly  linked  to  turbulent  stresses 
within  the  boundary-layer,  provide  a  more  significant  level  of  comparison  with  CFD  (Computational  Fluid 
Dynamics)  results  than  do,  eg,  surface  pressures.  Again,  however,  heat  transfer  and  especially  skin  friction 
are  difficult  to  measure  with  acceptable  accuracy  in  high-speed  interacting  flows. 

The  traditional  methods  for  the  measurement  of  skin  friction  in  high-speed  flows  have  been  thoroughly 
reviewed  by  Winter.4  Recent  developments  were  also  surveyed  by  Settles.3  Other  than  the  direct  measurement 
of  wall  shear  stress,  rtt,  by  a  balance,  the  available  techniques  infer  skin  friction  from  some  other  measured 
quantity  such  as  heat  transfer  or  pitot  pressure.  The  validity  of  such  inferential  methods  is  doubtful  in 
complex  interacting  flows.  Further,  floating-element  balances  themselves  have  serious  problems  in  such  flows, 
especially  due  to  pressure  gradients  and  poor  spatial  resolution.  Thus,  for  the  class  of  flows  in  question, 
reliable  measurement  techniques  have  been  essentially  nonexistent.  A3  a  result,  there  are  no  known  skin 
friction  data  in  shock/boundary  layer  interactions  of  accepatable  accuracy  for  code  validation  purposes. 

A  recent  development,  the  Laser  Interferometer  Skin  Friction  (LISF)  meter,  promises  to  resolve  this 
problem.  The  LISF  meter  was  invented  by  Tanner  and  Blows,0  and  was  subsequently  refined  by  Tanner, 
Monson  and  Higuchi,8  Monson,  Driver  and  Szodruch,2  Monson, 10,11  Westphal,  Bachalo,  and  Houser,12 
and  Kim  and  Settles.13  It  interferometrically  senses  the  time  rate  of  thinning  of  an  oil  film  on  a  polished 
surface  subjected  to  aerodynamic  shear.  In  two-dimensional  (2-D)  flows  without  pressure  gradients,  oil 
lubrication  theory  then  gives  rw  directly  without  any  reference  to  the  properties  of  the  overlying  boundary 
layer.  While  some  corrections  are  required  in  pressure-gradient  and  shear-gradient  flows,  the  instrument 
nonetheless  delivers  essentially  a  direct  measurement  of  rw.  The  Penn  State  version  of  this  instrument  and 
its  associated  data  acquisition  and  reduction  procedures13  is  believed  to  be  the  only  compressible- flow  LISF 
meter  currently  in  existence. 

The  applicability,  repeatability,  and  accuracy  of  the  LISF  technique  for  use  in  compressible  flows  have 
been  determined  by  Kim  and  Settles13  through  an  experimental  calibration  over  a  range  of  Mach  numbers.  It 
was  observed  that  surface-wave  phenomena  on  the  oil  film  limit  the  technique  in  high-r^  flows,  but  this  does 
not  preclude  accurate  measurements  in  the  supersonic  Mach  number  range.  Monson10  reached  a  maximum 
level  of  rw  =  120  N/m2  with  the  LISF  meter  in  high-Reynolds  number  supersonic  flow.  Here,  we  extend 
that  limit  to  almost  600  N/m2. 

The  present  paper  reports  an  experimental  study  in  which  the  LISF  meter  is  used  to  measure  tw 
distributions  in  two  interactions  of  fin-generated,  swept  shock  waves  with  turbulent  boundary- layers.  The 
basic  research  configuration  is  an  unswept,  sharp-leading-edge  fin  of  variable  angle,  a,  mounted  on  a  flat  plate. 
This  configuration  generates  “building-block”  interactions  representative  of  those  occurring  on  practical  high¬ 
speed  vehicles.  Two  examples  are  fin-fuselage  junctions  and  corner  flows  in  supersonic  engine  inlets. 

The  goal  of  this  study  is  to  apply  the  LISF  meter  to  obtain  reliable  rw  data  in  three-dimensional  (3-D) 
swept  shock  wave/turbulent  boundary-layer  interactions.  Such  data  are  sought  in  order  to  shed  light  on 
the  physical  mechanisms  of  such  interactions,  and  to  provide  a  benchmark  for  CFD  code  validation.  In 
particular,  the  data  are  expected  to  be  useful  in  judging  the  effectiveness  of  current  turbulence  models. 


Experimental  Methods 

Wind  Tunnel  and  Test  Conditions 

The  Penn  State  Supersonic  Wind  Tunnel,14  used  in  this  study,  is  a  blowdown  facility  with  a  nominal 
Mach  number  range  of  1.5  to  4.  The  variable  Mach  number  capability  is  achieved  by  way  of  an  asymmetric 
sliding-block  nozzle.  The  te3t  section  of  the  wind  tunnel  is  150  mm  (6  in.)  wide,  165  mm  (6.5  in.)  high,  and 
610  mm  (24  in.)  long.  Calibration  experiments  have  revealed  good  flow  quality  across  the  available  Macn 
number  range.  Wind  tunnel  run  durations  for  the  present  study  ranged  from  20  to  30  seconds. 

The  test  Mach  number  was  fixed  at  3.03 ±0.017  in  this  study.  The  stagnation  conditions  of  the  flow  were 
p0  =  0.827  MPa  (120  psia)  ±1.9%  and  Ta  =  293.5  K  ±1.4%,  yielding  a  freestream  unit  Reynolds  number  of 
6.19  x  107/m  ±2%. 

The  interaction  test  surface  was  a  flat  plate,  500  mm  (19.5  in.)  long,  which  spanned  the  test  section. 
The  plate  was  fitted  with  both  surface-pressure  taps  and  surface  thermocouples.  An  unswept  fin  shock- 
generator  was  positioned  with  its  sharp  leading-edge  216  mm  (8.5  in.)  from  the  plate  leading  edge  and 
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26.2  mm  (1.03  in.)  from  the  tunnel  sidewall.  The  Sn  was  100  mm  (4  in.)  high,  127  mm  (5  in.)  long  and 
6.35  mm  (0.25  in.)  thick.  The  fin  height  of  about  305  ensured  that  the  interaction  was  a  semi-infinite  one.1 
A  sketch  of  the  fin /plate  test  configuration  is  shown  in  Figure  1. 

The  fin  angle-of-attack  was  varied  by  a  pneumatic  fin-injection  mechanism  mounted  through  the  tunnel 
sidewalL  The  bottom  of  the  fin  had  a  rubber  3eal  which  ensured  that  no  leakage  under  the  fin  occurred 
during  the  experiments.  Two  fin  angle-of-attack  values,  a  =  10°  and  16°,  were  tested.  The  fin  angles  were 
determined  to  0.1°  accuracy  using  a  machinist’s  protractor.  -  ' 

The  flat  plate  contained  96  surface-pressure  taps  arranged  in  concentric  circu’ar  arcs  about  the  fin 
leading-edge  position.  This  arrangement  was  chosen  to  take  advantage  of  the  quasiconical  interaction  sym¬ 
metry  which  has  been  fully  described  elsewhere.1  The  pressure-tap  arc  radii  were  multiples  of  25.4  mm  (1 
in.  ).  LISF  measurements  were  also  made  along  circular  arcs  centered  about  the  fin  leading-edge.  However, 
these  were  chosen  midway  between  the  pressure-tap  arcs  in  order  to  avoid  interference.  For  the  or  =  10° 
case,  the  arc  radius  R  for  LISF  measurements  was  114.3  mm  (4.5  in.  ).  R  was  reduced  to  88.9  mm  (3.5  in.  ) 
for  a  =  16°  in  order  to  keep  the  measurement  area  away  from  possible  interference  due  to  the  fin  shock-wave 
interaction  with  the  wind  tunnel  sidewall  boundary-layer.  Finally,  some  additional  LISF  data  were  taken 
along  radii  from  the  conical-flow  virtual  origin.  These  radii  lay  at  0  =  24°  for  the  a  =  10°  case  and  /3  =  38° 
for  the  a  =  16°  case. 

The  undisturbed  flat-plate  boundary-layer  has  been  surveyed  and  documented14  for  test  conditions  close 
to  those  of  the  present  study  (same  Mx ,  but  Ree  11%  lower  than  present).  Pitot-pressure  surveys  along  the 
flat-plate  centerline  and  38  mm  (1.5  in.)  to  each  side  showed  that  the  test  boundarv-layer  is  two-dimensional 
turbuient,  and  in  equilibrium  in  the  sense  that  it  satisfies  the  combined  law-of-the-wall/Iaw-of-the-wake 
with  constant  IT  within  the  accepted  equilibrium  range.  The  test  boundary-layer  is  also  near-adiabatic 
{T.„/Tvv  =  1.06  ±  1.1%).  The  mean  velocity  profile  of  the  undisturbed  boundary-layer  at  x  =  28  cm 
(11  in.)  from  the  flat  plate  leading  edge  is  shown  in  Figure  2,  along  with  the  corresponding  Sun-Childs15 
compressible  wall-wake  curvefit.  The  boundary-layer  parameters  extracted  from  this  curvefit  are:  8  =  4.4 
mm,  Re^  =  13408,  IT  =  0.63,  and  C/  =  0.0015.  For  the  higher  Res  =  14915  of  the  present  study,  cj  in  the 
incoming  boundary-layer  is  predicted  to  be  0.0014  by  the  Van  Driest  II  theory. 

The  LISF  Meter  in  High-Speed  S-D  Flow 

The  underlying  theory  of  the  LISF  technique  and  examples  of  its  application,  mainly  in  incompressible 
flows,  are  thoroughly  documented  in  References  6-12.  The  adaptation  of  the  LISF  meter  for  use  in  compress¬ 
ible  flows  has  further  been  discussed  in  detail  by  Kim  and  Settles.13  The  major  limitation  of  the  technique 
is  due  to  surface-wave  phenomena  on  the  oil  film  in  high-r„,  flows,  about  which  more  will  be  said  below. 

Other  than  in  a  supersonic  delta-wing  demonstration  experiment  by  Monson,  et  al,9  the  LISF  meter  has 
not  been  used  in  compressible  3-D  flows  prior  to  the  present  study.  One  of  our  goals  is  co  firmly  establish 
the  applicability  of  the  instrument  to  such  flows. 

It  is  generally  necessary  in  3-D  LISF  measurements  to  know  the  direction  of  the  rm  vector  a  priori, 
since  an  accurate  knowledge  of  the  distance,  As,  from  the  oil-film  leading  edge  to  laser-beam  measuring 
spot  is  required.  Alternatively,  when  the  dual-beam  LISF  adaptation  is  used,8-12  As  need  not  be  measured. 
However,  in  that  case,  either  the  dual  beam  spots  must  still  be  aligned  with  the  local  rw  direction,  or  else 
two  components  of  rra  must  be  measured  with  considerable  extra  effort. 

Local  rw  directions  may  be  determined  directly  from  surface-flow  visualization  results,  assuming  that 
these  are  obtained  in  a  quantitative  manner.  The  kerosene-lampblack-adhesive  tape  technique16  is  particu¬ 
larly  suitable  for  this,  in  that  it  yields  undistorted  full-scale  surface-streak  patterns.  Angular  measurements 
of  local  tw  directions  are  possible  with  ±0.5°  routine  accuracy.  Such  patterns  were  obtained  for  the  a  =  10° 
and  16°  test  cases  of  the  present  study,  and  are  displayed  in  Figure  3. 

During  LISF  tests,  full-scale  transparent  overlays  of  these  surface-flow  patterns  were  used  to  identify  the 
surface-flow  direction  corresponding  to  a  given  position  of  the  laser-beam  measuring  spot.  After  applying 
the  oil  film  to  the  flat  plate,  an  optical  cathetometer  (least  count:  0.025  mm)  aligned  normal  to  the  surface- 
flow  direction  was  used  to  measure  As  with  a  repeatability  of  ±0.6%.  The  leading  edge  of  the  oil  film  was 
highlighted  for  this  purpose  by  reflected  light  from  a  spotlight  positioned  outside  the  wind  tunnel.  Still, 
uncertainty  of  the  exact  oil-film  leading-edge  position  was  the  major  contributor  to  uncertainty  in  As.  With 
practice,  As  was  measured  with  an  overall  accuracy  of  ±2%. 

With  this  technique  to  obtain  accurate  As  measurements  along  the  local  surface-flow  direction,  we  were 
able  to  use  the  single-beam  LISF  meter  arrangement  shown  in  Figure  1.  The  dual-beam  method,  originally 
developed  from  the  single-beam  method  by  Monson  and  colleagues,  was,  at  the  time,  a  solution  to  the 
problem  of  inaccuracy  in  measuring  As.  However,  for  3-D  flows  the  dual-beam  approach  appears  to  us  to 
involve  added  complexity  and  mechanical  alignment  difficulties  which  far  outweigh  its  advantages.  With  the 
simpler  single-beam  method,  there  is  no  basic  difference  in  the  LISF  technique  for  2-D  and  3-D  flows. 
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Arrangement  of  LISF  Meter  and  Data  Acquisition  System 

Figure  1  illustrates  the  components  of  the  single-beam  LISF  meter  used  in  this  study.  The  beam  from  a 
5-milliwatt  linearly-polarized  helium-neon  laser  (1)  first  passes  through  a  50%  neutral-density  filter  (2)  and 
an  iris  diaphragm  (3).  The  beam  is  then  directed  downward  by  a  folding  mirror  (4),  and  focused  by  a  lens 
(5),  passing  through  the  plexiglas  ceiling  window  (6)  of  the  wind  tunnel  to  form  a  spot  of  about  450  pm. 
diameter  at  an  appropriate  point  (7)  on  the  flat  plate.  During  expejynents,  a  thin  film  of  u  =  500  centistokes 
Dow-Corning  “200”  silicone  oil  applied  to  a  local  region  (7)  of  the-flat  plate  is  sheared  by  the  rw  distribution 
of  the  shock/ boundary-layer  interaction  due  to  the  fin  (8).  The  incident  laser  beam  is  reflected  by  both 
the  surface  of  the  oil  film  and  the  polished  plate  beneath  it.  This  produces  a  reflected,  two-component, 
interfering  laser  beam  directed  out  of  the  wind  tunnel  through  the  ceiling  v,  indow.  (The  beam  angles  of 
incidence  and  reflection  are  kept  to  about  1°.)  The  reflected  beam  is  intercepted  by  two  first-surface  aiming 
mirrors,  (9)  and  (10),  whence  it  is  directed  through  a  6328A  filter  (11)  onto  a  photodiode  (12). 

The  photodiode  senses  a  time-dependent  light  intensity  due  to  the  interference  of  light  reflected  from 
the  oil  film  and  test  surface.  The  photodiode  output  (13)  is  raised  to  a  level  of  8  volts  by  an  internal 
operational  amplifier,  low-pass  filtered  with  a  10  Hz  cutoff  to  remove  high-frequency  noise,  and  recorded. 
Seven  other  channels  of  data  are  recorded  simultaneously:  pQ,  T0,  two  pw  channels,  and  three  Tw  channels. 
The  9  data  channels  are  simultaneously  sampled  by  a  multiplexer,  digitized  by  a  microcomputer-controlled 
A/D  converter,  and  stored  in  the  computer  memory.  For  present  purposes  a  20Hz  data  rate  was  initially 
chosen,  yielding  at  least  20  points  per  interference  fringe.  However,  the  requirement  for  greater  accuracy  in 
high-shear  regions  ox  ihe  flow  led  to  an  increase  of  the  data  rate  to  50Hz.  The  photodiode  output  was  also 
monitored  on  a  stripchart  recorder  during  the  experiments. 

Data  Reduction 

A  detailed  account  of  LISF  data  reduction  procedures  is  given  by  Kim  and  Settles.13  Briefly,  the  limited 
fringe  count  obtainable  in  high-fu,  flows  calls  for  somewhat  different  data  reduction  than  that  reported  by 
other  investigators  for  incompressible  flows.  Given  only  2-8  fringes,  it  becomes  necessary  to  determine  tw 
based  on  the  entire  signal  within  the  usable  fringe  record.  Merely  determining  the  time  interval  corresponding 
to  3  arbitrary  fringe  peaks,  as  is  done  in  incompressible  flow,  does  not  yield  sufficient  accuracy  here. 

Thus,  as  described  in  Reference  13,  it  is  necessary  to  define  that  part  of  the  fringe  record  which  is 
usable  and  consistent  with  oil  lubrication  theory,  ie,  which  conforms  to  a  regular  “envelope.”  The  initial  and 
final  segments  of  the  fringe  record  which  do  not  so  conform  are  discarded,  as  are  occasional  distorted  fringe 
records  due  to  dust  particles  on  the  oil  film,  etc.  Usable  fringe  records  are  smoothed  by  repeated  application 
of  an  adjacent- point  averaging  algorithm  until  all  noise  on  the  signal  is  removed. 

The  LISF  data  reduction  equation,8  derived  from  the  Navier-Stokes  equation  and  assuming  constant 
wall  shear  stress,  is: 

2npi/cos(r)  [As  .  . 

=  A  [w]  (  } 

where  N'  =  effective  fringe  number,  t'  =  effective  oil  flow  time,  n  =  oil  refractive  index,  p  =  oil  density,  A  = 
laser  wavelength,  and  r  =  oil  refraction  angle.  For  present  purposes  this  equation  requires  a  correction  for 
the  effect  of  variable  Tw  on  t/,  as  described  in  Reference  13.  The  calculation  of  N'  and  t'  is  also  carried  out 
per  Reference  13,  eqns.  (6)  and  (7). 

These  results  aye  next  compared  with  oil  lubrication  theory  in  the  iV'-t'  plane,  as  shown  in  Figure  4  for 
a  case  where  rw  =  586  N/m2.  A  least-squares  curvefit  of  the  data  to  the  lubrication  theory  is  performed, 
followed  by  a  computation  of  the  reduced  x 2  deviation  of  the  data  from  the  fitted  curve.  Then,  the  first 
fringe  peak  (te,  the  first  data  point)  is  eliminated  and  the  fit  is  repeated.  This  process  is  continued  up  to 
A N  =  —1.5,  rather  than  -2  as  quoted  in  Reference  13,  considering  the  small  number  of  usable  fringes  in 
the  present  study  (due  to  high  rw)  compared  to  our  previous  work.  The  value  of  N't'  corresponding  to  the 
curvefit  with  the  minimum  reduced  \2  is  used  in  the  calculation  of  rw  in  Equation  1. 

Before  rw  is  found,  however,  two  further  corrections8,9  are  required  in  principle,  since  the  oil  film  in 
the  present  experiments  is  subject  to  both  pressure-  and  shear-gradients.  The  pressure-gradient  correction, 
per  Monson  and  Higuchi8,  was  evaluated  through  interpolation  of  pw  measured  at  the  radial  tap  rows  on 
either  side  of  the  LISF  measurement  arc  on  the  flat  plate.  These  pressure  distributions,  obtained  in  an 
earlier  study14,  are  plotted  in  Figure  5  as  functions  of  the  angle  /?.  Note  that  it  was  deliberately  chosen  to 
take  LISF  data  in  the  region  of  quasiconical  interaction  symmetry  well  downstream  of  the  fin  leading-edge. 
Thus,  as  shown  in  Figure  5,  pw  vs.  is  essentially  independent  of  R.  The  pressure-gradient  correction  for 
the  present  experiments  amounted  at  most  to  ±0.1%,  which  is  negligible.  This  fact  is  significant  in  that 
pressure-gradient  errors  are  controlling  factors  in  most  other  methods  of  skin  friction  measurement. 4,5 


Per  Monson,  Driver,  and  Szodruch,9  the  shear-gradient  correction  is  quite  simple  for  the  single-beam 
LISF  technique.  The  correction  is  required  if  there  is  a  significant  change  in  tw  from  the  oil-film  leading-edge 
to  the  laser-beam  measurement  spot.  This  correction  has  the  form: 


1  drw 

=  rw  +  -——As 
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where  t'w  =  the  corrected  value  of  wall  shear  stress.  For  small  As  one  may  assume  that  drw/ds  =  a  constant 
locally.  Thus  the  correction  b  accomplished  by  a  shift  of  the  measurement  point  along  the  surface  streamline 
direction.  However,  for  the  present  experiments  As  ce  2  mm,  so  the  required  shift  is  only  ~  0.5  mm,  which  is 
negligible.  In  general,  a  shear-gradient  correction  may  always  be  avoided  in  single-beam  LISF  measurements 
if  As  is  kept  small. 


Error  Analysis 

The  uncertainties  in  the  flow  quantities  and  other  parameters  connected  with  the  LISF  measurements 
have  been  cited  in  previous  sections.  In  all  cases  these  uncertainties  are  given  as  the  mean  value  ±  one 
standard  deviation.  Each  measurement  point  within  the  two  swept  interactions  studied  here  was  repeated  7 
times  in  separate  wind  tunnel  runs.  Chauvenet’s  criterion  was  then  applied  once  to  each  7-point  ensemble 
to  discard,  if  necessary,  a  single  bad  point  from  the  ensemble.  The  error  bars  in  the  plots  of  the  following 
sections  indicate  the  repeatability  of  che  data  as  ±  1  standard  deviation  about  the  mean.  In  general,  this 
repeatability  error  varied  from  as  little  as  ±2%  upstream  of  the  interaction  to  as  much  as  ±10%  at  the 
highest  rw  in  the  a  =  10°  case.  In  the  subsequent  a.  =  16°  experiments,  a  higher  data  acquisition  rate  and 
experience  combined  to  limit  the  repeatability  band  to  ±6%  or  less. 

Of  course,  there  exists  no  calibration  standard  with  which  to  judge  the  absolute  accuracy  of  tw  measure¬ 
ments  in  shock/boundary-layer  interactions.  However,  the  accuracy  of  LISF  measurements  in  compressible 
flat-plate  turbulent  boundary-layers  has  been  assessed  in  previous  calibration  experiments.9'13  By  coming 
to  terms  with  the  limits  imposed  by  high  shear  levels  on  the  usable  fringe  count,  the  error  band  has  been 
reduced  from  ±10%9  to  ±6% 13  at  the  Mach  number  of  the  present  study.  As  discussed  below,  further 
improvements  may  be  possible. 


Results  and  Discussion 

The  measured  skin  friction  distribution  for  a  =  10°,  R  =  114.3  mm  (4.5  in.)  is  shown  in  Figure  6.  As  in 
Figure  5,  the  »•/  data  are  plotted  vs.  the  angle  P  measured  from  the  fin  leading  edge,  since  the  interaction  is 
quasiconical.  The  flat-plate  c/  level  ahead  of  the  interaction  is  0.0011,  which  is  20%  below  the  value  predicted 
by  the  Van  Driest  II  theory.  At  the  upstream-influence  line  of  the  interaction  {p  =  39°),  c/  initially  rises, 
then  falls  again  in  the  region  between  the  separation  line  (p  =  32°)  and  the  angle  of  the  “inviscid”  fin  shock 
wave  {P  =  27.4°).  From  this  location  aft  to  P  =  15°  a  steep  rise  in  c/  occurs,  reaching  maximum  values  of 
0.0025  or  230%  of  the  incoming  level.  Finally,  c/  decreases  slightly  as  the  fin  (a  =  P  =  10°)  is  approached. 

Our  physical  interpretation  of  the  measurements  shown  in  Figure  6  is  as  follows.  The  initial  rise  of  c/  at 
the  upstream  influence  line  appears  to  be  due  to  the  onset  of  the  swept  pressure  gradient,  which  shears  the 
bottom  of  the  incoming  boundary-layer  in  the  spanwise  direction.  Low-momentum  fluid  is  thus  removed  from 
the  boundary-layer,  such  that  the  velocity-gradient  normal  to  the  surface  (and  thus  c /)  is  increased.  The 
onset  of  the  swept  interaction  may  therefore  be  imagined  as  a  natural  boundary-layer  control  mechanism. 
However,  separation  nonetheless  occurs  downstream,  where  the  convergence  of  upstream  and  downstream 
fluid  at  the  wall  reduces  the  normal  velocity-gradient  and  Cf .  The  steep  rise  to  maximum  c/  values  near  the 
fin  is  explained  in  terms  of  recent  structural  studies  of  separated  fin  interaction  flowfields.17  The  observed 
“A-shock”  structure  causes  the  impingement  of  a  supersonic  jet  on  the  flat  plate  between  the  fin  and  the 
line  of  flow  reattachment  upstream  of  the  fin.  This  phenomenon  is  known  to  produce  very  high  rw,  pw,  and 
heat  transfer  levels.18  From  this  jet  impingement  region  the  fluid  near  the  wall  moves  spanwise  and  forward 
to  eventually  meet  the  separation  line.  Thus  the  phenomenology  of  the  region  15°  <  P  <  32°  in  Figure  7  is 
a  rapid  decrease  of  Cf  from  its  maximum  level  as  the  separation  line  is  approached  from  behind. 

The  measured  3kin  friction  distribution  for  a  =  16°,  R  =  88.9  mm  (3.5  in.)  is  shown  in  Figure  7.  This 
stronger  interaction  exhibits  the  same  qualitative  c/  features  as  the  10°  case.  However,  the  maximum  cj 
level  is  now  0.0041,  or  370%  higher  than  the  incoming  boundary-layer  c/.  As  the  interaction  strength  grows, 
the  jet-impingement  phenomenon  described  above  also  grows  stronger  and  produces  higher  local  wall  shear, 
here  centered  about  P  ~  24°.  The  skin  friction  drops  rapidly  and  almost  symmetrically  on  either  side  of  this 
impingement  region. 


5 


IP  1  null  II,  » 1.1  r  I."  r  »J'  HM  HT  «.■■  ■  SJ'  Vi  WH'  -v  ■■.■»  vm  urn  >r»  »  -  ■  -  - 


The  Initial  rise  in  c/  at  the  upstream-influence  line,  prominent  in  Figure  6,  is  almost  unnoticeable  in 
Figure  7.  Here  the  upstream  influence  line  lies  at  fi  =  49°  and  the  separation  line  occurs  at  j)  =  45°,  ie,  very 
close  together.  Previous  studies1’19  have  shown  that  the  angular  increment  between  upstream  influence  and 
separation  decreases  with  increasing  interaction  strength,  eventually  becoming  quite  small.  Thus  the  initial 
rise  in  c;  at  the  interaction  onset  is  judged  to  be  primarily  a  weak-interaction  phenomenon. 

Between  the  inviscid  shock  location  (/3  =  33.3°)  and  the  separation  line,  c/  decays  less  rapidly  than  in 
its  steep  fall  from  the  jet-impingement  location  at  @  os  24°.  Withtf^thls  region  at  f)  =  37°  lies  a  surface-flow 
feature  commonly  referred  to  as  “secondary  separation,”  which  is  quite  weak  in  the  present  case  (see  Figure 
3b).  No  indication  of  its  presence  is  noticeable  in  the  c/  distribution. 

Throughout  these  experiments  the  performance  of  the  LISF  meter  was  quite  consistent  with  previous 
experience  in  simpler  flows.13  The  primary  noticeable  difference  in  performance  was  due  to  high  interaction 
shear  levels  which  preduced  fewer  usable  fringes  than  before.  However,  for  laser-beam  spot  locations  directly 
beneath  the  inviscid  fin  shock,  noisy  fringe  signals  were  obtained.  This  is  thought  to  result  from  laser-beam 
refraction  by  the  shock  when  the  two  are  essentially  coplanar.  Even  so,  it  was  possible  to  filter  out  the  noise 
and  still  obtain  a  Cj  measurement  beneath  the  shock. 

Comparison  with  Navier-Stokes  Computations 

Computational  simulations  have  been  carried  out  by  Horstman  and  Knight20  for  swept,  fin-generated 
interactions  with  essentially  the  same  boundary  conditions  as  in  the  present  experiments.  These  investigators 
provided  the  detailed  results  of  their  computations21'22  for  comparison  with  our  measured  c/  distributions.  In 
all  cases  the  compressible,  Reynolds-averaged  Navier-Stokes  equations  were  solved,  though  different  numeri¬ 
cal  algorithms,  grids,  and  turbulence  models  were  used.  Horstman’s  solutions  employed  the  Jones-Launder, 
2-eauation  k-e  turbulence  model,  while  Knight  used  the  Baldwin-Lomax  algebraic  eddy-viscosity  turbulence 
model.  A  more  detailed  description  of  these  computations  is  not  given  here  due  to  space  limitations,  but  is 
thoroughly  documented  elsewhere  (eg,  Reference  20). 

A  comparison  of  the  Horstman  and  Knight  CFD  results  with  the  present  a  =  10°  cj  data  is  shown 
in  Figure  6.  (Note  that  an  interpolation  of  the  computed  mesh  points  of  both  computations  was  done  by 
the  present  authors  in  order  to  present  all  results  at  R  =  114.3  mm  (4.5  in.)  from  the  fin  leading-edge.) 
Overall,  both  computations  reasonably  predict  the  high  measured  c/  level  around  /?  =  15°,  as  well  as  its 
decay  with  increasing  /?.  However,  the  initial  region  of  the  interaction,  from  the  upstream  influence  to  the 
inviscid  shock,  is  not  well  predicted  in  either  case.  Knight’s  result  shows  an  initial  decrease  in  C/,  whereas 
an  increase  was  measured.  Horstman’s  solution  does  predict  the  increase,  but  places  it  some  5°  aft  of  its 
measured  location.  Horstman21  recognizes  that  his  computational  simulations  of  such  interactions  tend  to 
underpredict  the  upstream  influence  for  reasons  unknown  at  present.  Knight’s  solution,  on  the  other  hand, 
matches  the  experimental  upstream  influence  angle  but  shows  the  opposite  trend  in  c/.  The  discrepancies 
in  the  incoming  boundary-layer  eg  level  in  Figure  6  are  due,  at  least  in  part,  to  somewhat  different  Reg 
boundary  conditions  among  the  3  cases  compared. 

A  comparison  of  the  Horstman’s  CFD  result  with  the  present  a  =  16°  cj  data  is  shown  in  Figure  7.  (No 
computation  for  this  fin  angle  was  carried  out  by  Knight.)  Again,  the  upstream  influence  is  underpredicted 
and  the  initial  cj  rise  is  overpredicted,  compared  with  the  measurements.  However,  Horstman’s  solution  is 
in  good  agreement  with  our  results  aft  of  the  inviscid  shock  location  at  /?  =  33.3°. 

Finally,  the  additional  LISF  data  taken  along  radii  from  the  quasiconical  virtual  origin  at  /3  =  24°  for 
the  a  =  10°  case  and  0  =  38°  for  the  a  =  16°  case  are  shown  in  Figures  8a  and  8b,  respectively.  In  Figure 
8a,  interpolated  values  of  both  the  Horstman  and  Knight  computations  are  also  shown.  The  abscissa  of 
Figures  8a  and  8b  is  given  as  the  streamwise  distance  from  the  virtual  conical  origin  of  the  interaction.1  This 
was  determined  from  surface-streak  patterns  in  the  experiment  and  from  the  computational  analogs  of  such 
patterns.  Presenting  the  results  in  these  coordinates  reveals  the  behavior  of  c/  in  the  inception  zone  of  the 
interaction  and  in  the  quasiconical  region  downstream. 

Figure  8a  reveals  that  cj  increases  with  distance  from  the  virtual  origin,  eventually  asymptoting  to  a 
constant  level.  The  measurements  and  both  computations  agree  on  this  point,  which  is  entirely  consistent 
with  the  quasiconical  interaction  model  proposed  earlier.1  Figure  8b  shows  the  same  phenomenon  in  the 
a  =  16°  interaction.  Horstman’s  computation  is  not  shown  here  because  its  underprediction  of  the  upstream 
influence  confuses  the  comparison  at  the  chosen  conical  ray. 

Implications  for  LISF  Measurements  in  High-Shear  Flows 

As  described  elsewhere, 9’13,23  surface-waves  on  the  oil  film  impose  the  most  serious  limitation  on  the 
LISF  technique.  By  disrupting  the  laser-beam  coherence,  these  waves  prevent  usable  interference  fringes 
from  being  obtained  until  the  oil  film  becomes  quite  thin.  This  effectively  limits  the  number  of  usable  fringes 
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one  can  obtain,  which  directly  affects  the  measurement  accuracy.  The  surface-wave  problem  i3  exacerbated 
as  tw  becomes  large. 

In  the  present  study  we  have  encountered  tw  levels  ranging  from  115  to  586  N/m2.  The  number  of 
usable  fringes  we  obtained  over  this  range  is  plotted  against  in  Figure  9.  This  Figure  indicates  that  an 
asymptotic  condition  of  roughly  2  usable  fringes  may  persist  even  at  extremely  high  shear  levels.  Actually, 
the  curvefit  to  the  data  shown  in  Figure  9  yields  the  expression: 
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but  extrapolation  of  this  fit  tc  higher  rw  values  is  not  warranted. 

The  upshot  of  Figure  9  is  the  following:  If  sufficient  accuracy  in  c/  can  be  had  from  only  2  interference 
fringes,  then  the  LISF  technique  can  still  be  used  despite  arbitrarily  high  shear  levels.  For  this  reason  we 
have  concentrated  our  efforts  here  and  in  Reference  13  on  extracting  the  maximum  useful  information  from 
an  extremely  limited  fringe  count.  There  appears  to  be  no  simple  way  to  increase  N,  although  switching  to 
a  green  helium-neon  laser  (A  =  5230A)  would  produce  a  slight  advantage.  The  fringe  count  could  be  at  least 
doubled  by  the  use  of  a  UV-laser,  but  this  approach  involves  several  serious  drawbacks. 

Conclusions 

Experimental  measurements  have  been  made  of  wall  shear  -tresses  in  swept  interactions  of  a  planar 
shock  wave  generated  by  a  fin  at  angle  of  attack,  and  the  two-dimensional  turbulent  boundary  layer  on  a 
flat  plate  at  Mach  3.03.  These  measurements  were  made  using  the  Laser  Interferometer  Skin  Friction  (LISF) 
meter.  The  significant  conclusions  of  this  study  are  as  follows: 

1.  The  LISF  technique  is  shown  to  be  practical  for  skin  friction  measurements  in  3-D  shock  wave/boundary- 
layer  interactions. 

2.  With  proper  data  acquisition  and  reduction  procedures,  the  LISF  technique  is  useful  even  at  extremely- 
high  wall  shear  stress  levels. 

3.  In  contrast  with  the  case  for  other  skin  friction  methods,  the  LISF  corrections  for  pressure-  and  shear- 
gradients  are  negligible  in  the  present  experiment. 

4.  Skin  friction  is  found  to  rise  at  the  onset  of  a  weak  swept  interaction,  and  to  reach  high  levels  near  the 
fin  for  both  weak  and  strong  interactions. 

5.  The  present  skin  friction  data,  obtained  in  7-point  ensembles  at  each  measuring  station,  has  an  overall 
repeatability  of  ±6%  or  better. 

6.  The  quasiconical  interaction  model  proposed  earlier  is  supported  by  the  present  results,  and 

7.  Navier-Stokes  computations  by  two  other  investigators  simulate  some  features  of  the  measured  cj  dis¬ 
tributions,  but  fail  to  predict  other  features  properly. 
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Fig.  7  -  cj  Distribution  for  a  =  16°  with  CFD  Com-  Fig.  8a  -  c{  Distribution  Along  0  =  24°  Conical  Ray 
parison.  for  a  =  10°. 
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Fig.  9  -  Number  of  Usable  Fringes  vs.  Shear  Level.  Fig.  8b  -  c/  Distribution  Along  =  38°  Conical  Ray 
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Abstract 

The  development  and  implementation  of  a  control  al¬ 
gorithm  for  the  microprocessor-based  stagnation  pressure 
control  system  of  the  Penn  State  Supersonic  Wind  Tunnel 
Facility  is  reported.  The  gas  dynamics  and  the  control- 
valve  characteristics  of  this  blowdown-type  facility  are  non- 
linearlv  related.  A  mathematical  model  was  developed  for 
the  open-loop  system  characteristics  and  was  linearized  for 
the  controller  design.  A  single-input,  single-output  PI  con¬ 
troller  was  chosen  for  this  task  herause  of  it-'  simplicity  and 
availability.  The  resulting  performance  of  the  supersonic 
wind  tunnel  was  found  to  be  quite  good,  with  stagnation 
pressure  variations  typically  held  to  within  1  to  2  percent. 

Nomenclature 

A  =  test  section  area,  m3 

.-1,  =  nozzle  throat  area,  m3 

C3  =  control  valve  ‘gas  sizing  coefficient'’ 

r„  =  specific  heat  at  constant  volume.  J/kg-K 

Cj,  i-  ;;  e-jflc  heat  at  constant  pressure.  J/kg-K 

hP  —  specific  entmdpy  r>t  air  through  valve,  ,1/kg 

h.  =  specific  enthalpy  of  air  through  rozzle,  J/kg 

Kp  =  controller  proportional  gain 

A",  =  controller  integral  time 

M  —  test  section  Mach  number 

m ,  =  storage  tank  air  mass,  kg  - 

m„  '  mass  flow  rate  through  valve,  kg/sec 

m„  =  stilling  chamber  air  mass,  kg 

m,  —  mass  flow  rate  through  nozzle  throat,  kg/sec 

P,l  —  desired  stagnation  pressure.  Pascal 

P,  =  storage  tank  air  pressure,  Pascal 

P„  =  stagnation  pressure.  Pascal 

R  —  gas  constant,  J/kg-K 

T|  =  storage  tank  air  temperature,  K 

T0  =  stagnation  temperature,  K 

L’i  =  storage  tank  air  internal  energy,  Joule 

U„  =  stilling  chamber  air  internal  energy,  Joule 

V’i  =  storage  tank  volume,  m3 

V0  =  stilling  chamber  volume,  m3 

v„  =  velocity  of  air  through  valve,  m/sec 

e.  =  velocity  of  air  through  nozzle  throat,  m/sec 

O  =  valve  opening  position,  deg 

p,  =  storage  tank  air  density,  kg/m3 

=  stilling  chamber  air  density,  kg/m3 
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Introduction 

This  paper  deals  with  the  solution  of  the  stagnation 
pressure  control  problem  in  the  Penn  State  Supersonic  Wind 
Tunnel  Facility.  A  sketch  of  the  facility  is  shown  in  Fig  1. 
This  facility  is  a  blowdown-type  wind  tunnel  with  an  asym¬ 
metric  sliding-block  nozzle  which  yields  a  variable  Mach 
number  capability  over  the  range  of  Mach  1.5  to  4.0.  Dry 
air  is  stored  in  a  56.6  m3  tank  at  up  to  2.07  MPa.  Testing 
requirements  call  for  regulated  stagnation  pressures  in  the 
range  of  0.34  MPa  to  1.72  MPa  for  typical  test  durations  of 
15  to  45  seconds. 

The  wind  tunnel  air  supply  system  is  equipped  with  a 
20.3  cm  diameter,  pneumatically-actuated  ball  valve.  If  this 
valve  is  held  at  a  fixed  position,  the  wind  tunnel  stagnation 
pressure  drops  during  testing  because  the  pressure  in  the 
storage  tank  is  continuously  decreasing.  Pressure  regula¬ 
tion  by  manually  adjusting  the  valve  opening  during  a  test 
is  possible  but  inaccurate.  An  obsolete  electromechanical 
pressure  controller  has  been  used  in  the  past,  but  with  un¬ 
satisfactory  performance  for  current  testing  requirements. 

For  example,  at  a  given  Mach  number,  it  is  sometimes 
required  to  maximize  the  test  duration  by  running  the  tun¬ 
nel  at  the  lowest  possible  stagnation  pressure  which  will 
maintain  t”’personic  flow.  Similarly,  it  is  often  necessary  to 
obtain  different  Reynolds  numbers  at  a  given  Mach  num¬ 
ber  by  way  of  preset  levels  "f  the  stagnation  pressure,  in 
all  cases,  a  constant  stagnate.,  pressure  during  testing  is  a 
basic  requirement. 

The  purpose  of  the  present  work  is  to  design  and  im¬ 
plement  a  controller  which  senses  the  stagnation  pressure 
and  adjusts  the  valve  automatically,  thus  maintaining  the 
stagnation  pressure  constant  regardless  of  test  Mach  num¬ 
ber. 

Control  systems  for  incompressible  fluids  have  been  dis¬ 
cussed  extensively  in  the  literature1,2.  However,  compara¬ 
tively  little  study  has  been  made  of  the  control  problems 
of  highly-compressible  flows.  The  major  difficulty  here  lies 
in  the  nonlinearity  of  both  the  gas  dynamics  and  the  vaive- 
nozzle  characteristics  in  the  case  of  a  supersonic  wind  tun 
nel.  Furthermore,  in  this  case  a  wide  variety  of  operating 
conditions  must  be  handled  as  described  previously.  Fi¬ 
lially,  a  highly-sophisticated  control  system  is  impractical 
from  considerations  of  economy  and  reliability  of  wind  tun¬ 
nel  operations. 

A  satisfactory  solution  to  this  control  problenrhas  been 
acliieved  through  a  single-loop  Proportional -plus- Integra l 
(PI)  controller  with  constant  parameter  settings.  The  first 
step  in  designing  this  control  system  is  to  develop  a  math¬ 
ematical  model  for  the  supersonic  wind  tunnel  flow,  as  de¬ 
scribed  below.  An  evaluation  of  various  control  algorithms 
is  then  presented,  along  with  the  rationale  for  the  selection 
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Fig.  1  The  Penn  State  Supersonic  Wind  Tunnel 
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of  the  most  suitable  one  for  present  purposes.  The  behavior 
of  the  control  system  has  been  simulated  by  computer  be¬ 
fore  its  actual  implementation.  Finally,  its  performance  has 
been  verified  bv  wind  tunnel  experiments,  whose  results  are 
compared  with  the  computer  simulation. 

Mathematical  Model 

A  prerequisite  for  the  control  system  design  is  the  de¬ 
velopment  of  a  mathematical  model  of  the  controlled  pro¬ 
cess,  which  is  the  flo>v  in  the  supersonic  wind  tunnel.  Two 
control  volumes  are  selected  for  this  analysis:  the  air  stor- 
nee  tank  and  the  wind  tunnel  stilling  chamber.  It  is  here 
assumed  that  all  thermodynamic  processes  are  isentropic 
during  a  wind  tunnel  test  (no  shock  waves,  frictional  ln.,ses. 
or  heat  transfer),  and  that  the  change  of  potential  energy 
of  the  gas  is  negligible.  In  reality,  it  is  recognized  that  only 
the  last  assumption  is  strictly  correct. 

Model  Equations  for  the  Storage  Tank 

Conservation  of  Mass  During  a  test,  it  is  assumed 
that  the  pressure,  temperature  and  density  distribution  are 
uniform  over  the  entire  control  volume.  Hence,  the  rate  of 
decrease  of  mass  in  the  air  tank  is  equal  to  the  rate  of  mass 
efflux  rhrough  the  control  valve: 

,,dPt 

v,7t  ~  ~mv 


-dr  =  -v,mv  (1) 

The  valve  characteristics  are  described  in  literature3'1 
by  its  manufacturer.  The  mass  flow  rate  at  different  valve 
positions  is  given  by: 

0.2958  x  10'”  „  (  fAP  \ 

m -  =  — It, —  Pt  sin  (2-71  V  T  )  .  (2) 


dU,  .  1  .  , 

-  -mvhv  -  -mvvv 


The  left-hand  side  of  this  equation  can  be  written  as: 

dU,  d  r_  .  ,,  :  d  \  Cv  Pt  V,  c„  V,  dP, 
dt  ~  dtLTn,  '~v^‘‘  ~  dt[  R  \~  R  dt 

Assuming  the  airflow  to  behave  both  calorically  and 
thermally  as  a  perfect  gas,  the  right-hand  side  becomes: 

rhvhv  +  -rhvi\2  -  mth,  =  mv  cf  T, 

Thus,  Equation  (3’)  can  be  written  as: 

C.  v,  dr, 

R  dt  ~~mv  °rTt 
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Model  Equations  for  the  Stilling  Chamber 

Conservation  of  Mass  Air  flows  into  the  stilling 
chamber  from  the  valve  connected  to  the  air  storage  tank, 
and  flows  out  through  the  converging-diverging  nozzle  to 
the  test  _ection.  Therefore,  the  rate  of  mass  influx  into 
the  stilling  chamber  control  volume  minus  the  rate  of  mass 
efflux  yields  the  net  rate  of  mass  buildup  in  the  stilling 
chamber. 


—  =  (4, 

The  inass  flow  rate  through  the  nozzle  is  a  function  of  stag¬ 
nation  conditions  and  nozzle  throat  area5,6: 

0.G847  P„  A. 


where  AP  is  the  pressure  difference  across  the  valve,  here 
assumed  to  be  given  by  AP  =  Pt  —  P„. 

Conservation  of  Energy  By  assuming  that  the  en¬ 
ergy  loss  through  the  control  valve  is  negligibly  small,  the 
internal  <*nergy  change  in  the  air  tank  is  set  equal  to  the 
enthalpy  and  kinetic  energy  flux  through  the  valve, 


For  a  specific  test  Mach  number  M  and  test  section  area  A, 
the  nozzle  throat  area  is  given  by: 

,  1.728  A  M 

A.  =  - 5-  (6) 

(1  +  0.2  A/2)3 


Conservation  of  Energy  The  energy  entering  the 
stilling  chamber  from  the  valve  with  mass  m„  minus  the 
energy  exiting  through  the  nozzle  with  mass  m.  is  equal  to 
the  rate  of  increase  of  internal  energy  in  the  stilling  chamber, 

dL  s  .  ,  1  2  ,  1  .  2  /— ( 

— r-  =  m„h„  +  -mvvv  -  m.h.  -  -m.u.  (i  ) 
dt  2  2 

The  left-hand  side  can  be  written  as: 

d£ 

dt 
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Rewriting  the  right-hand  side  of  Eqn.  (?')  in  terms  of  rp\ 
mvhv  +  =  m,  ht  =  mv  cp  Tt 


m.h.  4-  -rh.v.2  =  m.  h„  =  m.  cp  T0 

Therefore,  the  statement  of  energy  conservation  for  the  still¬ 
ing  chamber  becomes: 

c„  Vo  dP,  .  . 

~  “  W  y  Cp  i  (  “  Ml  a  Cp  1  (J 


R  dt 
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Mathematical  Model  for  the  Wind  Tunnel 

From  the  preceding  discussion,  these  expressions  were 
obtained  which  describe  the  characteristics  of  the  airflow  in 
the  .  up-  i sonic  wind  tunnel: 

dpt 


dt  vt  m" 


dP « 
dt 
dPo 
dt 


pp  R  \  -7- 

-Z-rr  mv  T, 

Cv  ►</ 


dP i  {  cpR\  ,  ,  ■  -r  i 

!T  =  lzvJ 


(i) 

(3) 

(4) 
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where  the  mass  flow  rate  through  the  valve  and  the  nozzle 
are  given  by  Eqns.  2,  5  and  6. 

Expressing  temperature  as  a  function  of  pressure  and 
density  via  the  perfect-gas  equation  of  state,  the  above  equa¬ 
tions  become  four  first-order  nonlinear  ordinary  differential 
equations  in  the  four  state  variables  P, ,  pt,  P0  and  p0.  The 
inputs  to  this  system  are:  test  Mach  number  M,  which 
specifies  the  nozzle  geometry,  and  valve  position  0,  which 
determines  the  valve  behavior.  The  output  of  the  system 
is  the  stagnation  pressure  P0,  The  system  of  equations  was 
solved  for  P ( \f.  0 )  by  way  of  a  standard  finite-difference 
computer  routine. 

In  order  to  judge  the  accuracy  of  this  mathematical 
model,  experimental  and  simulated  results  for  stagnation 
pressure  vs.  time  at  Mach  4  with  the  control  valve  fully 
open  are  compared  in  Fig.  2.  The  initial  air  pressure  of 
the  storage  tank  in  this  comparison  is  1.8  MPa  (260psia). 
The  startup  process  of  the  wind  tunnel  involves  about  a 
2-second  delay  which  is  not  simulated  by  the  model.  Other¬ 
wise.  however,  the  model  accurately  predicts  the  observed 
rate  of  decay  of  P„  with  time. 


HO6  Math  eo 


0  5  LO  15  20  25  30 

Time  (see) 


Fig.  2  Stagnation  Pressure  vs.  Time,  Valve  Fully  Open 

Controller  Design 

Control  Algorithm  Analysis 

Upon  establishing  the  mathematical  model  for  the  su¬ 
personic  wind  tunnel,  a  control  algorithm  based  on  this 
model  was  sought.  The  goals  of  the  control  action  of  this 
algorithm  are: 

(1)  To  obtain  a  constant  steady-state  stagnation  pressure 
with  as  little  transient  oscillation  as  po^rible, 

(2)  To  maintain  this  constant  pressure  as  long  as  possible, 
and 

(3)  To  meet  the  above  requirements  at  any  desired  com¬ 
bination  of  test  Mach  number  and  stagnation  pressure 
within  the  ranges  cited  previously. 

Several  candidate  control  algorithms  were  considered 
for  this  task,  including  optimal  control,  state  feedback  con¬ 
trol  and  single-loop  P1D  control. 

In  the  optimal  control  of  a  process,  a  control  decision  is 
made,  subject  to  certain  contraints.  that  will  minimize  some 
measure  of  deviation  from  the  ideal  performance7.  The  per¬ 
formance  index  of  the  system  evolves  from  a  consideration 
of  such  factors  as  stability,  efficiency,  economy,  reliability, 
etc.  Once  the  performance  index  is  defined,  the  best  pos¬ 
sible  system  is  that  which  minimizes  this  index,  subject  to 
the  contraints  imposed  upon  the  system. 

In  short,  the  present  wind  tunnel  control  problem  is  rel¬ 
atively  simple  in  that  only  accuracy  and  stability  are  mat¬ 
ters  of  prime  concern.  It  was  thus  judged  that  the  complex¬ 
ities  of  optimal  control  are  neither  needed  nor  desired  for 
present  purposes. 

The  state-space  control  method  is  based  on  the  de¬ 
scription  of  the  system  in  terms  of  n  first-order  differential 
equations8.  In  the  present  rase,  Eqns.  1,  3,  4,  and  7  are  the 
state-space  equations  of  the  dynamic  system,  and  P,,  P„,  p, 
and  p„  are  the  state  variables  of  this  system.  The  design 
of  a  full  state  feedback  controller  involves  the  mathematical 
description  of  the  desired  system  response,  followed  by  a  de¬ 
termination  of  the  feedback  gain  of  each  state  variable  that 
yields  this  response.  At  present,  however,  only  two  state 
variables,  P,  and  P„,  are  directly  measurable.  Real-time 
calculation  of  pt  and  p0  would  thus  be  required,  rendering 


this  approach  cumbersome  ami  potentially  expensive  to  im¬ 
plement.  Once  again,  a  simpler  approach  is  called  for,  eg, 
PID  control. 

Tlie  principle  of  the  PID  control  scheme  is  to  act  upon 
the  variable  under  control  through  a  proper  combination  of 
three  control  actions:  Proportional,  Integral,  and  Derivative 
control9.  Proportional  control  is  proportional  to  the  actu¬ 
ating  error  signal  (  the  difference  between  the  input  and  the 
feedback  signal)  itseif.  while  integral  and  derivative  control 
ate  proportional  to  the  integral  and  derivative  of  tiie  actu¬ 
ating  *»rror  signal,  respectively.  For  the  present  sii  person*'* 
wind  tunnel  application  of  PID  control,  the  input  to  the 
controller  is  the  desired  stagnation  pressure  P,{,  while  the 
feedback  signal  is  the  actual  stagnation  pressure  P3.  This 
feedback  control  process  is  represented  by  a  block  diagram 
shown  in  Fig.  3. 


Fig.  3  Block  Diagram  of  PID  Control  Process 


is  Tint  obviously  needed.  Thus,  for  simplicity,  we  consider 
only  a  proportional- plus- integral  (PI)  controller  at  present. 
The  possible  need  for  additional  derivative  control  may  be 
judged  from  the  expeiimental  results  to  be  presented  later. 
The  transfer  function  of  a  PI  controller  is  given  by: 


G(s)  =  —  =  Kp  (l  +  — ) 

E(3)  pV  Kia) 


(3) 


where  .1  stands  lor  the  time-derivative  in  the  Laplace  domain 
a::!  F's)  is  the  error  signal  between  the  reference  input 
Pi  (desired  stagnation  pressure)  and  the  system  output  F„ 
(the  actual  pressure  measured  in  the  stilling  chamber).  The 
n.ove  equation  can  be  rewritten  to  yield: 

0(s)  =  Kp  (l  +  [Pd  -  Po(3)j,  or 

tei.  =  kp .,  +  1^.  {pi-  /»,(,))  (9) 

Applying  the  inverse  Laplace  transform  in  Eqn.  9,  the  dif- 
'crential  relations! jp  He'we“n  'he  input  E{t)  =  Pj  —  P ,'•) 
and  output  0(t)  of  the  PI  controller  is: 


Tli  •  effects  of  proportional,  integral  and  derivative  con¬ 
trol  action  on  the  performance  are  briefly  described 

as  follows:  Proportional  control  is  the  algorithm  in  which 
the  change  in  the  control  signal  is  proportional  to  the  devi¬ 
ation  of  'he  system  output,  from  t he  deseed  value-0.  <5clo 
proportional  control  of  a  process  introduces  a  steady-state 
ertor.  or  oifset.  in  the  response  to  a  step  input  This  is  the 
cliief  disadvantage  of  a  purely-proportional  controller. 

The  integral  control  of  a  process  is  equivalent  to  adding 
an  integrator  to  an  open-loop  system11,  integral  control 
action  at  any  instant  is  proportional  to  the  area  under  the 
actuating  error  signal  curve  up  to  that  instant.  As  long 
a-s  there  is  an  error  signal,  the  integral  controller  produces 
an  output  to  reduce  this  error.  Therefore,  the  steady-state 
error  of  a  proportional  control  system  can  be  eliminated  if 
integral  control  is  added. 

Derivative  control  action  responds  to  the  rate  of  change 
of  the  actuating  error.  It  thus  “anticipates"  the  error  and 
initiates  an  early  corrective  action.  Where  required,  this 
tends  to  increase  the  transient  response  and  stability  of  the 
system. 

Overall.  PID  control  is  a  well-developed  method  for  a 
single-input  single-output  system.  Only  one  measurement 
of  a  system  output  is  needed  for  the  control  process.  The 
simplicity  and  reliability  of  PID  control  are  its  major  ad¬ 
vantages  over  other  methods.  Moreover,  it  is  usually  the 
least  expensive  choice  of  those  discussed  above.  These  ad¬ 
vantages  led  to  the  choice  of  the  PID  process  for  the  solution 
of  the  present  supersonic  wind  tunnel  control  problem. 

Controller  Design  Analysis 

The  transients  at  startup  and  shutdown  of  the  super¬ 
sonic  wind  tunnel  are  quite  brief.  Thus  the  major  concern  in 
controller  design  is  focused  on  the  steady-state  performance 
of  the  wmd  tunnel.  In  order  to  achieve  and  maintain  a 
preselected  value  of  stagnation  pressure,  both  proportional 
and  integral  control  actions  are  clearly  required.  How¬ 
ever,  the  transient-handling  capability  of  derivative  control 


(16(f) 
c it 


=  —  A', 


dPM)  I< 


dt 


+ 


A',- 
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where  the  term  dPp{t)/dt  is  known  from  the  relation  of 
energy  conservation  of  'he  stilling  chamber  (Eqn.  7). 

The  open-loon  system  ( ie  without  feedback  control)  is 
highly  nonlinear,  both  in  terms  of  the  gas  dynamics  and  the 
control  valve  characteristics.  Following  standard  control- 
theory  procedure,  the  system  was  linearized  at  selected 
Mach  numbers  and  stagnation  pressures  over  the  operat¬ 
ing  range  of  the  wind  tunnel  facility.  This  was  done  by  way 
of  Tayior-scries  expansions,  dropping  all  higher-order  terms. 

Next,  the  controller  parameters  for  different  operating 
points  were  determined  as  follows: 

(1)  A  desired  closed-loop  system  with  specified  pole  loca¬ 
tions  is  established  first. 

(2)  The  PI  control  action  (Eqn.  10)  is  added  to  the  lin¬ 
earized  model  for  different  operating  conditions. 

(3)  The  linearized  model  together  with  the  control  action 
is  fitted  to  the  desired  closed-loop  system  to  determine 
suitable  controller  parameters  Iip  and  K,.  A  proper 
combination  of  I(p  and  A',-  relocates  the  poles  of  the 
linearized  open-loop  system  to  the  desired  locations  de¬ 
termined  in  step  (1)  above. 

Values  of  Kp  and  A',  were  thus  found  for  each  desired 
operating  condition.  These  values  were  observed  to  vary  by 
no  more  than  about  20%  over  the  entire  wind  tunnel  Mach 
number  and  stagnation  pressure  range.  Thus,  as  a  first  ap¬ 
proximation,  constant  average  values  of  Kp  =  0.00014  and 
A ,  =  0.08  sec-1  were  assumed  for  all  operating  conditions. 
If  this  approach  later  proves  not  to  be  sufficiently  accu¬ 
rate.  the  controller  parameters  may  be  allowed  to  vary  as 
functions  of  the  wind  tunnel  operating  conditions,  using  the 
concept  of  gain  scheduling5. 

Closed-Loop  System  Mathematical  Model 

The  PI  control  action  (Eqn.  10)  together  with  the  or;g- 
innl  open-loop  system  mathematical  model  (Eqns.  1,  3,  4, 


and  7)  yield  a  set  of  expressions  for  the  complete  feedback 
control  system.  Combining  Eqns.  7  and  10,  the  control  ac- 
'ion  of  the  PI  controller  takes  the  following  form: 


dO  r.  (  cr  R 


{ m.Ta  -  mvT,)  -f  -j?-  [P*  -  Pa)  (11) 

^  I 


Thus  the  complete  mathematical  description  of  the  closed- 
loop  system  is  described  concisely  by  Eqns.  1,  3,  4,  and 
7,  representing  the  gas  dynamics  of  the  supersonic  wmd 
tunnel,  and  Eon.  11  representing  the  control  action  of  the 
PI  controller.  The  desired  stagnation  pressure  Pd  in  Eqn  11 
is  the  reference  input  of  this  closed-loop  system. 

Analysis  of  Closed-Loop  System  Simulation 

As  was  done  for  the  open-loop  wind  tunnel  equations 
described  previously,  the  above  mathematical  model  for  the 
closed- loop  system  has  been  solved  by  a  finite-difference 
computer  routine.  The  resulting  pressure-‘ime  curves  for 
both  the  storage  tank  and  stilling  chamber  for  test  section 
Mach  numbers  of  2.5  and  4.0  are  shown  in  Figs.  4a  and 
4  u  r*v?p»’Ctive:y.  The  valve  opening  position  vs.  Mme  is  albo 
shown  for  "he  Mach  2.5  case  in  Fig.  4c. 
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Fig.  4a  Controlled  PQ  vs.  Time  at  Mach  2.5 
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Fig.  4c  Valve  Opening  Position  vs.  Time  at  Mach  2.5 


Analysis  of  the  simulation  results  reveals  the  following: 

(1)  The  constant  average  controller  parameters  found 
above  (Kp  =  0.00014  and  K,  =  0.08  sec-1)  are  ef¬ 
fective  at  all  Mach  numbers  in  obtaining  a  response 
with  a  minimum  steady-state  error  and  overshoot  as 
well  as  a  minimum  settling  time.  These  controller  set¬ 
tings  were  also  found  to  be  capable  of  maintaining  the 
different  steady-state  stagnation  presstires  constant  for 
all  Mach  numbers  in  the  range  of  2.5  to  4.0.  Thus  the 
only  required  input  to  the  closed-loop  system  is  the 
desired  reference  pressure  Pj. 

(2)  Maximum-duration  testing  is  possible  by  presetting  Pj 
to  the  minimum  value  to  sustain  supersonic  flow  at 
a  given  Mach  number.  Similarly,  variable-Reynolds 
number  testing  is  possible  over  a  wide  P„  range  at  a 
given  Mach  number.  Of  course,  the  simulation  predicts 
shorter  test  durations  for  higher  P0  values.  However, 
both  these  modes  of  testing  are  well-simulated  using 
the  fixed  controller  parameters  found  above. 

(3)  At  a  given  Mach  number  and  a  given  reference  input, 
the  effects  of  adjusting  the  controller  parameters  on  the 
system  response  are  as  follows: 

•  If  the  proportional  gain  is  increased,  the  maximum 
pressure  overshoot  and  the  settling  time  are  in¬ 
creased.  But  at  the  same  time,  the  steady-state 
response  can  be  maintained  for  a  longer  time. 

•  If  the  proportional  gain  is  decreased,  oscillations 
of  P0  during  the  initial  transient  period  are  re¬ 
duced,  but  the  steady-state  test  time  is  also  re¬ 
duced.  This  is  because  a  higher  proportional  gain 
makes  the  system  response  more  oscillatory,  but  it 
also  responds  more  readily  to  the  dropping  tank 
pressure,  especially  near  the  end  of  the  test  when 
a  large  increase  in  valve  opening  is  called  for. 

•  Changing  the  integral  time  within  a  certain  range 
does  not  have  an  obvious  influence  on  the  closed- 
loop  system  response. 


Fig.  4b  Controlled  P,  vs.  Time  at  Mach  4.0 
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Having  developed  a  successful  theoretical  model  for  the 
supersonic  wind  tunnel  and  its  stagnation  pressure  con¬ 
troller.  we  now  consider  the  hardware  aspect  of  the  problem. 
A  suitable  commercially-available  microprocessor  PID  con¬ 
troller  (Moore  Products  Co.  MYCRO-35210)  was  obtained 
and  installed  in  the  console  of  the  Penn  State  Supersonic 
Wind  Tunnel.  While  not  designed  specifically  for  this  ap¬ 
plication,  this  controller  is  programmable  with  built-in  func¬ 
tion  blocks  which  may  be  -.cufigured  to  yield  a  wide  variety 
of  different  control  algorithms. 

The  microprocessor  controller  communicates  with  the 
wind  tunnel  via  built-in  A/D  and  D/A  converters.  The 
former  accepts  and  digitizes  an  amplified  signal  represent¬ 
ing  the  “process  variable  from  a  pressure  transducer  which 
senses  Pa  in  the  stilling  chamber.  The  latter  sends  a  4mA 
to  20  mA  analog  output  signal  to  a  “reverse"  (current-to- 
pressure)  transducer  which  modulates  the  air  pressure  to  the 
bellows- type  control  valve  actuator.  Single-input,  single¬ 
output  operation  is  thus  achieved,  with  all  controller  actions 
either  *aking  place  automatically  within  the  microprocessor 
controller  or  resulting  from  manual  input  to  the  controller 
by  the  wind  tunnel  operator. 

Prior  to  a  wind  tunnel  test,  the  operator  preselects  the 
desired  stagnation  pressure  P4  and  enters  a  corresponding 
“set  point"  on  the  display  of  the  controller.  The  “set  point" 
and  the  “process  variable"  provide  inputs  to  the  internal 
PID  controller  function  block  of  the  MYCRO-352,  which 
compares  these  to  find  the  actuating  error  signal,  then  cal¬ 
culates  and  produces  an  appropriate  output  signal  to  the 
control  valve. 

Other  operator  inputs  to  the  microprocessor  controller 
include  A’p,  A',,  and  the  derivative  gain  (here  set  to  zero). 
Manual  valve  control  is  also  available  to  the  operator  hv 
way  of  an  automatic/manual  switch.  An  auxilliary  on/off 
switch  in  the  controller  output  line  is  used  to  initiate  and 
terminate  wind  tunnel  runs.  For  safety  reasons,  the  control 
valve  actuator  is  set  up  such  that  either  a  loss  of  signal  from 
the  controller  or  a  loss  of  electrical  power  in  the  laboratory 
will  result  in  immediate  valve  closing. 

The  panel  display  of  the  microprocessor  controller  pro¬ 
vides  the  wind  tunnel  operator  with  information  including 
the  set  point,  process  variable,  and  vaive  opening  indices. 
A  bar- graph  display  allows  visual  comparison  of  the  actual 
vs.  preselected  stagnation  pressure  levels.  Various  alarm 
functions  are  available  for  programming  if  desired.  Finally, 
a  pulser  knob  is  provided  for  manual  input  to  the  controller 
by  the  operator.  Further  details  on  the  controller  and  its 
programming  may  be  found  in  Refs.  10  and  11. 

Experimental  Verification  and  Discussion 

The  computer  simulation  of  the  wind  tunnel  and  con¬ 
trol  system  has  been  compared  with  the  actual  operation  of 
the  Penn  State  Supersonic  Wind  Tunnel.  In  so  doing,  the 
simulation  served  as  an  important  guide  in  understanding 
the  relationship  between  control  action  and  the  correspond¬ 
ing  closed-loop  system  response,  and  in  allowing  the  trial- 
and-error  to  be  done  by  computer  instead  of  in  the  actual 
wind  tunnel. 

A  comparison  of  the  stagnation  pressure  vs.  time  of  the 
model  simulation  and  the  win  !  tunnel  test  data  is  shown  in 
Figs,  5a-c  for  Mach  2.5,  3.0,  and  3.5.  The  closed-loop  system 
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Fig.  5b  Comparison  of  Simulation  and  Data  at  Mach  3.0 


Fig.  5c  Comparison  of  Simulation  and  Data  at  Mach  3.5 


response  of  the  supersonic  wind  tunnel  is  very  similar  to  the 
results  predict  ’d  by  the  model  simulation.  In  every  case  a 
constant  controlled  P0  was  obtained  in  close  agreement  with 
the  model.  The  few  features  of  the  real  feedback  control 
system  not  predicted  by  the  simulation  are  discussed  below: 
( I !  Tile  same  constant  controller  parameters  determined  in 
the  simulation  were  used  in  t lie  experiments  and  were 
found  to  he  suitable  for  all  Mach  numbers  in  the  range 
of  2.,'j  to  3.5  at  different  stagnation  pressures.  However. 
Mach  numbers  below  2.3  were  not  attempted  in  the 
experiments  due  to  tunnel  blockage  caused  by  the  flat- 
plate  model  installation  in  the  test  section. 

(21  In  tiie  final  seconds  of  a  run.  the  valve  should  open  very 
rapidly  to  maintain  a  constant  P0.  In  reality,  however, 
the  vnive  actuator  cannot  respond  at  the  required  rate 
for  this  to  occur.  This  actuator  lag  is  not  accounted 
for  in  the  simulation.  Thus,  in  actual  testing,  when 
the  tank  pressure  is  only  20-30  psia  above  the  steady- 
state  stagnation  pressure,  both  pressures  will  decrease 
gradually. 

i  3'  For  'he  test  results  reported  here,  no  second-fhroa*  dif¬ 
fuser  was  installed  in  the  wind  tunnel  facility.  The 
minimum  running  stagnation  pressures  are  thus  rather 
high,  ranging  from  80  psia  at  Mach  2.5  to  220  psia  at 
Mach  4.0.  Only  quite  short  runs  were  possible  at  Mach 
4.0  in  this  case,  since  the  initial  difference  between  the 
tank  pressure  and  P„  was  at  most  60  psia.  The  recent 
installation  of  a  fixed  second  throat  iias  increased  this 
difference  to  100  psia  or  more,  providing  a  more  reason¬ 
able  Mach  4.0  testing  capability  of  at  least  30  seconds 
at  constant  P„.  No  change  in  controller  programming 
was  required  in  this  case. 

( 4  I  Tiie  filter  breakpoint  frequency  of  the  A/D  converter  in 
the  MYCRO-352  controller  is  10  Hz.  This  means  that 
the  controller  is  unaware  of  any  change  in  the  feed¬ 
back  signal  between  1/10  sec  sample  intervals.  Fince 
tiie  controller  output  is  a  step  function,  the  closed-loop 
system  response  is  adjusted  every  0.1  second,  giving  a 
slightly  oscillatory  driving  output  to  the  valve  actuator. 
It  is  unlikely  that  the  actuator  can  respond  at  this  fre¬ 
quency.  Small  oscillations  do  occur  in  the  measured  P„ 
of  the  wind  tunnel,  but  these  amount  at  most  to  only 
2  psia  for  P„  -  IGOpsia  (1.2  %),  which  is  considered 
acceptable. 

Summary 

Achieving  constant  stagnation  pressure  is  a  critical  con¬ 
cern  for  supersonic  blowdown  wind  tunnel  testing.  In  the 
present  study,  a  control  algorithm  has  been  developed  to 
acliieve  this.  The  control  algorithm  is  designed  such  that 
it  is  suitable  for  different  Mach  number  testing  and.  at  the 
same  time,  obtaining  the  maximum  test  time  for  different 
stagnation  pressures. 

A  mathematical  model  based  on  tiie  thermodynamic 
properties  of  air  was  established  for  the  supersonic  wind 
tunnel.  Comparison  of  the  system  response  of  the  model 
simulation  and  tiie  actual  wind  tunnel  test  data  was  made 
to  determine  the  applicability  of  the  model.  Once  a. suitable 
system  model  was  found,  it  was  then  used  to  design  the 
stagnation  pressure  controller. 


After  investigating  different  control  algorithms,  a 
single-input,  single-output  PI  controller  was  chosen  for  this 
task  because  of  its  simplicity  and  availability.  The  major 
problem  in  implementing  this  control  system  is  the  high  de¬ 
gree  of  nonlinearity  of  both  the  gas  dynamics  and  the  control 
valve  characteristics.  A  linearized  mathematical  model  was 
used  to  analxze  the  open-loop  system  characteristics  as  well 
as  for"the  controller  design. 

A  microprocessor -based  TID  controller  was  installed  in 
the  Penn  State  Supersonic  Wind  Tunnel  and  programmed 
according  to  the  computer  prediction.  The  resulting  perfor¬ 
mance  of  the  wind  tunnel  at  different  Mach  numbers  and 
stagnation  pressures  was  quite  good,  with  the  stagnation 
pressure  typically  being  held  constant  to  within  1%  to  2% 
over  a  wide  range  of  testing  conditions. 
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